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Preface 


The mining industry behaves a significant part in the current society. They provide 
coal and hard rock resources to guarantee the normal operation of the current society. 
To guarantee that coal mines and hard rock mines can operate smoothly, ground 
control techniques are compulsory. Among kinds of ground control techniques, rock 
bolts are widely accepted and applied in mine sites. They can improve the stability 
of the excavated rock masses. Consequently, the safety of underground excavations 
and operators can be guaranteed. 

In situ tests show that although rock bolts have been used for a long time, failure 
of the rock bolting system still occurs. Moreover, it is more common to encounter the 
bond failure of the bolt/grout interface. This leads to improper load transfer between 
rock bolts and surrounding rock masses. Furthermore, this induces safety concern 
for the normal operation of mine sites. 

Aiming at solving this issue, this book summarises the recent research work 
regarding by load transfer of rock bolts conducted by the authors. This book is 
divided into seven chapters. Within those seven chapters, the shear behaviour of the 
bolt/grout interface, the shear failure of the cement-based grout, the load transfer 
performance of rock bolts and the parameter analysis on the performance of rock 
bolts are studied. This book is beneficial for researchers and engineers to better 
understanding the load transfer performance of rock bolts. Moreover, it can help 
researchers and engineers to propose and develop new approaches to prevent failure 
of the rock bolting system. 

The publication of this book is financially supported by the National Natural 
Science Foundation of China (52034009, 52174093, 51904302). The authors also 
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would like to thank the support provided by the China University of Mining and Tech- 
nology (Beijing), and International Joint Research Laboratory of Henan Province 
for Underground Space Development and Disaster Prevention, Henan Polytechnic 
University, China. 


Beijing, China Jianhang Chen 
Yongliang Li 
Junwen Zhang 
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Chapter 1 A) 
Introduction E 


1.1 Background 


Mining engineering plays a significant role in the current society. It provides the 
energy resources including coal and hard rocks to guarantee the normal operation of 
the society. Therefore, not only developed countries but also developing countries 
still pay much attention to their mining industry. 

This is more obvious in China. For example, the statistical data shows that more 
than a half of the electricity still comes from the thermal power generation. Therefore, 
continuous production of coal resources is important in guaranteeing the normal 
power generation and convenient running of the society. 

To make sure that coal mines can operate stability, the ground control techniques 
are significant. Itis well known that after the rock masses around the mineral resources 
are excavated, the rock masses around the excavation tend to move towards the exca- 
vation direction. Therefore, to maintain the stability of the underground excavations, 
ground control techniques and instruments should be used. 

Among kinds of ground control techniques and instruments, rock bolts are widely 
accepted and used by mines. The rock bolt is a bar or a cable that is installed in the 
rock mass. Although rock bolts can be divided into frictional bolts or mechanical 
bolts, the laboratory and in situ tests show that the mechanical bolts have better 
performance. 

For mechanical rock bolts, either resin-based or cement-based grout can be used 
as the anchorage agent. These different grouts can be poured into the drilled borehole 
to bond the rock bolt with the surrounding rock mass. As for the in situ practice, the 
point-anchored method or the fully grouted method can be used. However, previous 
research shows that when the fully grouted method is used, there is a better load 
transfer between the rock bolt and the surrounding rock masses [1]. 

Nevertheless, when the fully grouted rock bolts are used, failure of the rock bolting 
system still occurs. A summary of the rock bolting failure modes shows that the 
rock bolting system can fail by different types, including the bond failure at the 
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bolt/grout interface, the bond failure at the grout/rock interface, shear failure of the 
grout column, shear failure of the rock mass and tensile rupture of the rock bolt [2]. 

However, it is more common to encounter the bond failure at the bolt/grout inter- 
face [3]. The reason is that the shear strength of the bolt/grout interface is much 
smaller than the tensile strength of the rock bolt. Moreover, in the rock bolting 
system, the bolt/grout interface has the relatively smaller contact area. This leads to 
the shear stress concentration at the bolt/grout interface. Consequently, bond failure 
of the bolt/grout interface is more likely to occur. 

To better understand the load transfer mechanism of fully grouted rock bolts and 
prevent bond failure of the bolt/grout interface, this book summarised the previous 
research work conducted by the authors regarding the anchorage performance of rock 
bolts. This book is divided into seven chapters. Within those seven chapters, the shear 
behaviour of the bolt/grout interface, the shear failure of the cement-based grout, the 
load transfer performance of rock bolts and the parameter analysis on the performance 
of rock bolts are studied. This book is beneficial for researchers and engineers to 
better understanding the load transfer performance of rock bolts. Moreover, it can 
help researchers and engineers to propose and develop new approaches to prevent 
failure of the rock bolting system. 


1.2 Book Outline 


This chapter aims at providing the general introduction for this book. The general 
background of rock bolting is illustrated. Then, the basic structure of this book is 
presented. 

Chapter 2 gives a literature review on the shear stress propagation mechanism 
in the rock reinforcement system. More attention is paid to the bolt/grout interface. 
The investigation approaches used in previous research are summarised. Moreover, 
the shear stress propagation process at the bolt/grout interface is elaborated. And a 
detailed discussion is added to analyse the shear stress propagation mechanism of 
the bolt/grout interface. 

Chapter 3 studies the shear behaviour of the cement-based grout. Two different 
water-cement (w/c) ratios are used: 0.4 and 0.5. The numerical compressive strength 
tests are conducted to calibrate the input parameters. Then, the direct shear test is 
simulated to study the shear behaviour of the cement-based grout. Experimental 
direct shear tests are used to confirm the accuracy of the numerical direct shear test. 
Moreover, the shear stress distribution in the cement-based grout is analysed. 

Chapter 4 presents an analytical modelling to investigate the shear behaviour of the 
bolt/grout interface. The bond-slip behaviour of the bolt/grout interface is simulated 
with a tri-linear model. An analytical model is developed based on the shear stress 
propagation mechanism at the bolt/grout interface. The pull-out performance of the 
rock bolts is divided into the elastic, elastic-softening, elastic-softening-debonding, 
softening-debonding and debonding stages. Credibility of the analytical model is 
validated with experimental pull-out tests. 
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Chapter 5 studies the influence of the confining medium on the load transfer 
performance of rock bolts. In this chapter, the tri-linear model is still used to depict 
the shear behaviour of the bolt/grout interface. Then, this tri-linear model is incor- 
porated into the rock bolting system. After that, based on this analytical model, the 
influence of the confining medium on the load transfer capacity of rock bolts is 
studied. Moreover, the critical influence diameter of the confining medium for the 
rock bolts is analysed. 

Chapter 6 conducts a comprehensive parameter study on the load transfer perfor- 
mance of rock bolts. The influence of the bolt diameter, elastic modulus of bolts, 
grouted length of bolts and bond slipping when bond strength reaches is studied. The 
maximum load transfer capacity of rock bolts is analysed and compared. Moreover, 
when the rock bolt reaches the maximum load transfer capacity, the shear stress 
distribution state of the bolt/grout interface is analysed. 

Chapter 7 proposes an analytical model for fully grouted rock bolts based on a 
nonlinear bond-slip model. In this chapter, an exponential nonlinear model is used to 
depict the shear behaviour of the bolt/grout interface. Then, the load transfer perfor- 
mance of rock bolt is successfully simulated. Experimental pull-out tests are used to 
confirm that this nonlinear model can be used to simulate the pull-out performance 
of rock bolts. Moreover, the load distribution along the axial direction of rock bolts 
is studied. 
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Chapter 2 A) 
Analysis on the Shear Stress Propagation s 
Mechanism in the Rock Reinforcement 

System 


2.1 Introduction 


In rock mechanics, failure of rock masses is commonly encountered [1]. This is 
because fractures distribute non-uniformly in rock masses [2, 3]. Moreover, experi- 
ment work proves that rock masses which are full of fractures have quite low strength 
[4]. 

This is more prominent when the rock masses are subjected to manual excavation 
[5]. For example, in civil engineering and mining engineering, excavation activi- 
ties are performed to create tunnels, chambers or roadways [6]. These tunnels and 
roadways will be later used in serving transportation and ventilation [7]. 

Attention is paid that this manual excavation disturbs the initial stability of rock 
masses [8]. Moreover, due to manual excavation, stress concentration occurs in rock 
masses [9, 10]. Consequently, fractures are likely to develop in rock masses [11-13]. 
This further weakens the rock mass strength. 

To guarantee the safety of rock mass excavation, rock reinforcement bolts are 
commonly used [14—16]. Experimental work proves that rock reinforcement bolts 
can combine the jointed rock mass. Moreover, once rock masses converge, shear 
deformation will occur in the grout column. Consequently, stress can be transferred 
in the rock reinforcement system [17]. 

During the loading process, the shear stress propagation plays a significant role 
[18]. Therefore, proper understanding the shear stress propagation mechanism in the 
rock reinforcement system is quite significant. 

This paper aims at revealing the shear stress propagation mechanism in the rock 
reinforcement system. To realise this purpose, a literature review was conducted. 
This study is beneficial to propose new reinforcement approaches to prevent rock 
mass failure. 
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2.2 Previous Study on the Shear Stress Propagation 


2.2.1 Investigation Approaches 


To study the shear stress propagation mechanism in the rock reinforcement 
system, previous researchers adopted various approaches. Generally, investigation 
approaches can be classified into three different types: experimental tests [19], 
analytical simulation [20, 21] and numerical simulation [12]. 

Experimental tests should be the most credible one. Specifically, strain gauges are 
adhered on the bolt/grout interface [22]. Then, this instrumented bolt is installed in 
a rock block. Grout is used to bond the bolt with the rock block [23, 24]. After full 
curing, the bolt is pulled out. During testing, tensile force distribution along the bolt 
can be measured. Then, Eq. (2.1) is adopted to simply calculate the shear stress at 
the bolt/grout interface [25]. 


oS Dp Ep(€1 — €2) 


4AL ep 


where t: shear stress at the bolt/grout interface; Dp: bolt diameter; Ep: elastic modulus 
of the bolt; e: tensile strain of the bolt; AL: spacing between two adjacent strain 
gauges. 

This approach has been used by a number of researchers [26]. However, a short- 
coming is that the attached strain gauges are likely to be stripped from the bolt/grout 
interface. To solve this issue, Chekired et al. [27] developed the tension measuring 
device. This device can be mounted on the bolt to measure the strain distribution. 
Additionally, Martin et al. [28] proposed replacing central wire in the cable bolt with 
a modified wire. Specifically, along this modified wire, strain gauges were attached. 

Besides, analytical simulation is an efficient approach to study the shear stress 
propagation mechanism [29]. Specifically, although the bolt may have a long length, 
infinitesimal method can be used to analyse the shear stress at the bolt/grout interface, 
as shown in Eq. (2.2) [30]: 


Dp do (x) 
T= —> 
4 dx 


(2.2) 


where do,(x): the increment of the tensile stress in the bolt. 

Then, it is assumed that the bolt/grout interface has the same mechanical property 
[31]. And a bond-slip equation can be used to depict the relationship between the 
shear stress at the bolt/grout interface and the slip [32]. Through incorporating the 
bond-slip equation into Eq. (2.2), analytical equations can be developed. 

This analytical approach was initially proposed by Farmer [33]. However, at that 
period, the debonding behaviour was not considered. Later, this analysis method was 
adopted by others [34]. Aydan et al. [35] made a significant improvement on analytical 
simulation. Their innovation is that the classic tri-linear equation was proposed to 
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describe the slip behaviour of the bolt/grout interface. Based on the tri-linear equation, 
the bolt/grout interface encountered the elastic, softening and debonding behaviour. 
Therefore, the bonding and debonding behaviour of the bolt/grout interface can be 
simulated. 

Moreover, the tri-linear equation may be simplified into the bi-linear equation 
to evaluate the shear stress propagation of the bolt/grout interface [36]. The main 
difference between the bi-linear equation and the tri-linear equation is that there is 
no linear softening behaviour in the bi-linear equation. 

In recent years, numerical simulation becomes more popular in rock reinforce- 
ment analysis. It is because numerical simulation is powerful in establishing and 
calculating complicated structures [37]. As for the rock reinforcement analysis, more 
research was focused on using the structure element, such as the cable or pile. There 
is a significant difference between the cable and pile. Specifically, the cable only 
considers the longitudinal performance of rock reinforcement bolts. In contrast, the 
pile can analyse both the longitudinal performance and the lateral performance of 
rock reinforcement bolts. 

The advantage is that numerical elements have already been created by the 
commercial company. Therefore, users can conveniently adopt the structure element 
to simulate different rock reinforcement cases. In contrast, the shortcoming is that the 
original constitutive equation may not truly reflect the proper behaviour of the rock 
reinforcement system. For example, in the structure element of cable, the bolt/grout 
interface is assumed to deform following an elastic perfectly plastic equation, as 
shown in Eq. (2.3) [38]. 


T = koô (ô < êp) 
(=a a 


where kg: shear stiffness of the grout column; ô: slip of the bolt/grout interface; 6): 
slip of the bolt/grout interface when the peak strength reaches; tp: peak strength of 
the bolt/grout interface. 

It neglects the post-failure behaviour of the bolt/grout interface. Therefore, it 
cannot truly simulate the loading performance of bolts without modification. 

Nevertheless, the commercial software usually reserves the secondary develop- 
ment interface. For example, for the structure element of cable, the Itasca Company 
prepares a number of FISH functions [39]. For the structure element of pile, the Itasca 
Company creates the TABLE function [40]. Therefore, users can use these FISH func- 
tions or TABLE function to modify the original rock reinforcement elements. For 
example, Fig. 2.1 shows the comparison between the original pile and the revised 
pile. Apparently, with the original pile, at the end of the shear stress propagation 
process, the shear stress at the full bolt/grout interface equalled the peak strength. 
This overestimated the loading capacity of bolts. In contrast, with the modified pile, 
there was always a non-uniform shear stress distribution at the bolt/grout interface. 
This was more consistent with the experimental test results. Consequently, it saw a 
wide application of the commercial numerical tools in rock reinforcement analysis 
[41]. 
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Fig. 2.1 Comparison between the structure element of pile in analysing the shear stress propagation 
process: a original pile; b revised pile 


2.2.2 Shear Stress Propagation Mechanism 


Based on previous research, it is accepted that the shear stress at the bolt/grout 
interface may have a uniform distribution if the anchor length was short enough 
[42]. Benmokrane et al. [43] indicated that when the anchor length is less than four 
times of the bolt diameter, the shear stress can be treated equal. Based on this concept, 
Eq. (2.4) was used to calculate the shear stress at the bolt/grout interface [44]: 


F 
y a 
x Dy L 


(2.4) 
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where F: pull-out force; Z: anchor length. 

Attention should be noted that Eq. (2.4) is valid when the anchor length is 
constant during the pull-out process. Nevertheless, it is more common to encounter 
the scenario where the anchor length decreases gradually. Then, Eq. (2.5) can be 
used to calculate the shear stress [44]. 


F 
pet ee (2.5) 
x Dy(L — up) 


where m: pull-out displacement. 

More importantly, rock reinforcement bolts usually have a long length [45]. In 
this case, after the bolt is subjected to tensile loading, non-uniform shear stress 
distribution occurs. For the laboratory monotonous loading condition, the bolt usually 
has two ends. One is embedded in the rock block, and it is called as the internal end. 
By contrast, the other one is left outside and it is called as the external end. Since 
there is a non-uniform shear stress distribution, shear stress propagation between two 
ends of the bolts generates. 

Some research indicated that the maximum shear stress under each loading level 
was likely to occur around the same position [46]. Moreover, that position was close 
to the borehole collar. 

In contrast, it is more common to see that during the initial load process, shear 
stress at the borehole collar increased gradually. With the loading increasing, the 
shear stress at the borehole collar increased to the peak strength. Then, it started 
dropping. More interestingly, the maximum shear stress moved towards the internal 
end direction. This shear stress propagation ended when the shear stress at the internal 
end of the bolt reached the peak strength. 

As a validation of this shear stress propagation concept, Rong et al. [47] conducted 
laboratory pull-out tests on bolts. Strain gauges were attached on the bolt/grout inter- 
face to record the tensile force distribution. Later, Ma et al. [44] analysed these exper- 
imental data. Equation (2.1) was used to calculate the shear stress at the bolt/grout 
interface. The analysis results showed that at the initial loading grade, the shear 
stress at the borehole collar increased gradually. However, after a certain loading 
level, the shear stress at the borehole collar reached the peak strength. With the 
loading level further increasing, the maximum shear stress propagated gradually 
towards the external end. This analysis result was consistent with the above shear 
stress propagation concept. Therefore, the experimental work and the corresponding 
data analysis proved the reliability of the shear stress propagation concept. 

Additionally, the analytical simulation and numerical simulation can better reflect 
the shear stress propagation process. Ren et al. [48] used the classic tri-linear equa- 
tion to depict the slip behaviour of the bolt/grout interface. They analysed the shear 
stress propagation process at the bolt/grout interface. The results showed that the 
maximum shear stress at the bolt/grout interface propagated from the external end to 
the internal end. Moreover, although each point at the bolt/grout interface experienced 
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the elastic, softening and debonding behaviour, the full bolt/grout interface under- 
went five different grades. They were the elasticity, elasticity-weakening, elasticity- 
weakening-debonding, weakening-debonding and debonding grades. Later, Blanco 
Martin et al. [49] indicated that when the tri-linear equation was used to depict the 
slip behaviour of the bolt/grout interface, the full bolt/grout interface may undergo 
the pure softening grade. 

It should be mentioned that when different bond-slip equations are used, the full 
bolt/grout interface may undergo different grades. For example, Chen et al. [50] 
indicated that when a bi-linear equation was used to depict the slip behaviour of the 
bolt/grout interface, the full bolt/grout interface only experienced three grades: the 
elastic grade, elastic-debonding grade and debonding grade. 

Although different bond-slip equations can be used, the shear stress propagation 
mechanism was consistent. Specifically, the maximum shear stress at the bolt/grout 
interface consistently propagated from the external end to the internal end, as shown 
in Fig. 2.2. 

This finding was also confirmed with numerical simulation. Nemcik et al. [38] 
modified the original structure element in Fast Lagrangian Analysis of Continua 
(FLAC) and simulated the shear stress propagation process at the bolt/grout interface. 
A nonlinear bond-slip equation was used. The results showed that each point at the 
bolt/grout interface obeyed the same nonlinear bond-slip equation. With the loading 
level increasing, the maximum shear stress at the bolt/grout interface propagated 
towards the internal end. This finding was consistent with the others [51]. 


2.3 Discussion 


The shear stress at the bolt/grout interface plays a significant role in determining the 
loading capacity of bolts [52]. Under the static loading condition, the loading capacity 
of bolts equals the sum of the shear force at the bolt/grout interface [53]. Since the 
shear stress at the bolt/grout interface can be calculated directly with the shear force 
at the bolt/grout interface, as shown in Eq. (2.6), there is a close relationship between 
the shear stress at the bolt/grout interface and the loading capacity of bolts [54]. 


a (2.6) 
T= A 


where F,: shear force at the bolt/grout interface; As: contact area between the bolt 
and the grout column. 

Therefore, it is valuable to understand the shear stress propagation mechanism 
of bolts. To realise this purpose, previous researchers used various approaches. It is 
believed that the experimental approach is more credible. This is because compared 
with the experimental approach, the analytical simulation and numerical simulation 
usually relied on a number of assumptions [55]. Whether those assumptions are 
reasonable is doubted. In analytical simulation, it is usually assumed that the bolt, 
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grout column and surrounding rock masses deform elastically. In fact, in experimental 
tests, failure of the rock masses may occur because of the radial dilation at the 
bolt/grout interface. Therefore, under this scenario, the analytical simulation which 
assumes that only slip failure occurs at the bolt/grout interface may not be trustable. 

Similarly, in numerical simulation, rock masses are simulated with different 
elastoplastic equations. Among kinds of equations, the Mohr—Coulomb equation 
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is more commonly used. This equation is relatively simpler, and its input parame- 
ters can be acquired directly from experimental tests. However, it cannot properly 
simulate the post-failure behaviour of rock masses. Specifically, after the peak, the 
strain softening behaviour and residual behaviour of rock masses cannot be properly 
simulated. In this case, the interaction between the numerical rock mass and bolts 
cannot be the same as the reality. This is also the reason why the numerical simulation 
work should be calibrated and compared with experimental results. 

Nevertheless, this is not to deny the significance of analytical simulation and 
numerical simulation. In fact, these two approaches are more powerful for researchers 
to understand the shear stress propagation mechanism of bolts. Therefore, it is 
suggested that combing the experimental tests, analytical simulation and numer- 
ical simulation is a better choice in studying the shear stress propagation mechanism 
of bolts. 


2.4 Conclusions 


This paper conducted a literature review on the shear stress propagation mechanism of 
rock reinforcement bolts. The previous investigation approaches were summarised. It 
is concluded that previous researchers usually used the experimental tests, analytical 
simulation and numerical simulation. Among those approaches, the experimental 
approach is more widely used. And the experimental test results are more likely to be 
accepted by others. By contrast, the analytical approach had more degree of freedom. 
With this approach, researchers can use different bond-slip equation to depict the slip 
behaviour of the bolt/grout interface. As for the numerical simulation, it is convenient 
for users since the original constitutive equation has already been created by the 
developer. Moreover, developers usually reserve the secondary development interface 
for users to modify the original constitutive equations. 

As for the shear stress propagation mechanism, it is more commonly agreed that 
during the pull-out process of bolts, the shear stress at the borehole collar firstly 
increased. With the loading level increasing, the shear stress at the borehole collar 
gradually reaches the peak strength. Then, with the further increasing of the loading 
level, the maximum shear stress starts propagating towards the internal end. This 
phenomenon was consistently observed in experimental tests, analytical simulation 
and numerical simulation. 
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Chapter 3 A) 
Numerical Simulation of the Shear E 
Behaviour of Cement Grout 


3.1 Introduction 


Cement is grey fine powder that is widely used in civil engineering and mining engi- 
neering [1, 2]. In civil engineering, as an important ingredient, cement is always 
mixed with water, sand and rocks to form concrete, which is the commonly used 
construction or building material [3, 4]. In mining engineering, cement is usually 
mixed with water to form the cement grout, which is then poured into drilled bore- 
holes to bond the reinforcing tendon, such as rock bolts and cable bolts, with the 
surrounding rock masses [5]. Moreover, the cement can be used as the grouting 
material to be poured into rock mass to improve its mechanical behaviour. 

The mechanical properties of the cement grout have a significant effect in deciding 
the performance of the rock reinforcement systems [6]. Consequently, a vast of 
research has been conducted to study the mechanical behaviour of the cement grout 
[7]. Those studies can be basically divided into studying the axial behaviour of the 
cement grout and shear behaviour of the cement grout. 

To study the axial behaviour of the cement grout, researchers mainly used the 
Uniaxial Compressive Strength (UCS) tests. Feldman and Beaudoin [8] indicated 
that the strength and elastic modulus of the Portland cement grout were largely 
influenced by the porosity. Domone and Thurairatnam [9] found that with the w/c 
ratio decreasing, the strength of the cement grout increased. Hyett et al. [10] evaluated 
the influence of curing time on the grout strength and found that the grout UCS 
increased nonlinearly with the hydration time ranging from 1 to 7 days. Boumiz 
et al. [11] adopted the ultrasonic and calorimetric methods on cement grout. They 
found that the larger the hydration degree, the higher the elastic modulus of the 
cement grout. Li et al. [12] studied the influence of nanoparticles on the performance 
of cement grout. They proved that adding nanoparticles was beneficial for improving 
the grout strength. Mirza et al. [13] revealed that the grout strength increased within 
the curing period, independent of the cement type. 
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The above research played an important role in revealing the mechanical perfor- 
mance of cement grout. However, experimental tests proved that in the fully grouted 
rock reinforcing system, failure of the system usually occurred at the bolt/grout inter- 
face [14]. This was resulted by the low shear strength of the interface [15]. The shear 
behaviour of the cement grout has a significant role in determining the shear strength 
of the bolt/grout interface [16]. Therefore, better understanding the shear behaviour 
of the cement grout was beneficial for preventing failure of the bolt/grout interface. 

To study the shear behaviour of cement grout, researchers commonly used tri- 
axial tests and direct shear tests. For tri-axial tests, Hyett et al. [17] studied the 
performance of Portland cement grout with the confining pressure ranging from 5 
to 15 MPa. The results showed that the larger the confining pressure, the higher 
the strength of the cement grout. Moreover, the cohesion and the internal friction 
angle were successfully acquired. Simms and Grabinsky [18] performed tri-axial 
tests on cement material under the confined condition. The results showed that in the 
shearing process, the cementitious material had an apparent tendency to dilate. Sakai 
et al. [19] studied the performance of the cement paste under the tri-axial condition. 
The results proved that increasing the confining pressure changed the sample failure 
mode. Specifically, when the confining pressure was 30 MPa, oblique shear failure 
plane occurred in the sample. However, once the confining pressure was much larger, 
horizontal failure plane occurred in the sample. Aziz et al. [20] conducted tri-axial 
tests on two types of cement grout and used Mohr—Coulomb (MC) model to fit the 
shear strength envelope. A vast of tri-axial tests were conducted by Chen and Xu 
[21] to study the performance of the cement grout under the high confining pressure 
condition. It was found that with the confining pressure increasing, the performance 
of the cement grout changed from brittle to ductile. Moreover, this phenomenon 
was valid independent of the w/c ratio. Those tri-axial testing research improved the 
understanding in the shear behaviour of cement grout. However, in tri-axial tests, the 
shear failure plane occurred dependent on the properties of the grout and the loading 
status [22]. This was not a true reflection of the shear failure of the bolt/grout interface 
in the rock reinforcing system, where the shear failure plane was a pre-defined plane 
[23, 24]. 

Therefore, the direct shear test was a better approach to study the shear behaviour 
of the cement grout for the rock reinforcement scenario. The direct shear test is a 
traditional method to evaluate the shear strength of a material. Normally, a shear 
box is prepared. Then, the tested sample is cut to the desired geometry with the 
required dimension. After that, the sample was set in the shear box and installed in 
the apparatus. Operators can apply normal stress on the sample and then apply shear 
velocity along the shear direction. This test is continued until the shear strength was 
obtained. In this way, the shear strength of the material can be acquired. Moosavi and 
Bawden [25] designed a shear box to conduct direct shear tests on Portland cement 
grout. The results showed that the direct shear test was likely to provide a smaller 
shear strength of the cement grout, when compared with tri-axial tests. Lu and Wang 
[26] performed direct shear tests on cementitious samples with two different w/c 
ratios: 0.3 and 0.4. They found that when the w/c ratio was 0.3, shear strength of 
the cement-based sample increased apparently with normal stress. However, when 
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the w/c ratio was 0.4, shear strength of the cement-based sample did not increase 
much with the normal stress. Assaad et al. [27] proposed using the direct shear test to 
study the yield stress of cement material. The results revealed that the material shear 
strength was largely dependent on the normal stress, following the MC rule. Chen 
et al. [28] compared the shear performance of a Portland cement grout under the 
constant normal load and constant normal stiffness condition. It was found that the 
boundary condition might influence the performance of the cement grout. Li et al. 
[29] adopted the direct shear tests to investigate the influence of particle size on the 
shear strength of cement grout. They found that through adding superfine cement, 
the early-age shear strength parameters of the cement grout can be largely enhanced. 

That research provided new insights and knowledge in revealing the shear perfor- 
mance of cement grout. However, most research was conducted with experimental 
tests. Compared with that, much less work has been performed via numerical simu- 
lation. Therefore, this study aims at studying the shear behaviour of the cement grout 
with numerical simulation. 

First, a numerical calculating programme was selected to conduct this research. 
After that, numerical direct shear tests were performed to evaluate the shear behaviour 
of the cement grout. Then, the numerical calculating results were validated with 
experimental results. Last, the shear stress distribution on the shear failure plane was 
investigated. 


3.2 Selection of the Software 


In this study, the three-dimensional software FLAC3D was used. This programme 
was developed by the Itasca Consulting Group. 

The reason to use this programme is that the Itasca Company developed the 
interface element in it. This interface element is composed of triangular elements 
and corresponding nodes. Moreover, the interface element is defined with the shear 
coupling constitutive model. In this case, once shear displacement occurs in the 
interface element, shear stress at the interface accumulates and can be calculated. 
Relied on this function, the shear behaviour of a material can be studied. 

The authors of this study previously used the two-dimensional software FLAC 
to simulate the shear failure process of the bolt/grout interface [30]. However, the 
drawback was that the problem was simplified to be two dimensional. Compared 
with FLAC, the advantage of FLAC3D was that it can simulate objects in three 
dimensions. Therefore, it can better reflect the mechanical behaviour of objects. 
Then, in this study, the FLAC3D was used to simulate the shearing process of the 
cement grout. 
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3.3 Numerical Simulation Process 


A vast of direct shear tests were performed by Moosavi and Bawden [25]. Those tests 
were used in this study and simulated with FLAC3D. In their study, two different 
w/c ratios were used, namely 0.4 and 0.5. Then, in the following section, numerical 
tests on those two cement grouts were performed. 


3.3.1 Simulation of the Cement Grout with a Higher Strength 


3.3.1.1 Numerical UCS Test and Validation 


Before the direct shear test was conducted, UCS tests were performed to acquire the 
basic mechanical properties of the grout, including the UCS, elastic modulus and the 
Poisson’s ratio. Detailed information about the UCS test process and results were 
reported by Moosavi [24]. Experimental results showed that for the cement grout 
with a w/c ratio of 0.4, the UCS, elastic modulus and Poisson’s ratio were 50.6 MPa, 
12.66 GPa and 0.25. More specifically, in the UCS test, the cement grout was filled 
in PVC tubes with the diameter of 50 mm and the length of 150 mm. 

In numerical simulation in FLAC3D, the w/c ratio cannot be directly simulated. As 
a compromise, the UCS of the numerical mesh was selected as the index to compare 
with the experimental result. In FLAC3D, anumerical UCS test was created (Fig. 3.1). 

In the numerical UCS test, a cylindrical mesh was established. At the top and 
bottom of the cylinder, two square steel plates were established to simulate the 
loading plates. The whole mesh was composed of 16,032 zones and 16,531 grid 
points. For the numerical plate, there were 32 zones and 100 grid points. For the 
numerical sample, there were 16,000 zones and 16,441 grid points. 

A strain softening model was used for the cylindrical sample. It was developed 
based on the MC model. Compared with the MC model, the strain softening model 
had additional input parameters: the cohesion table and the internal friction angle 
table. With those two parameters, the cohesion and internal friction angle of the 
materials can decrease with the increasing plastic shear strain. Therefore, this model 
can simulate the post-peak behaviour of materials. 

The reason to use this strain softening model was that this model can simulate the 
post-peak behaviour of materials. Therefore, it is better in reflecting the mechanical 
behaviour of materials. 

An elastic model was used for the square plates. The reason to use the elastic 
model to simulate the square plate was that in the experimental compressive test, the 
cement-based sample was compressed by the steel plate, since the strength of the 
steel plate was much higher than the cement-based sample. Failure only occurred 
in the cement-based sample. The steel plate consistently kept the elastic state. The 
input parameters for the cylindrical sample and the square plates are tabulated in 
Table 3.1. 
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Fig. 3.1 Geometry of the UCS test created in numerical simulation 


Table 3.1 Input parameters 


for the sample and the plate Croup Input parameter yaus 

Sample Bulk modulus (GPa) 8.44 
Shear modulus (GPa) 5.06 
Cohesion (MPa) 19.3 
Internal friction angle (°) 21 
Tension (MPa) 2 

Plate Bulk modulus (GPa) 133.33 
Shear modulus (GPa) 80 


The velocity of the plates was fixed. Moreover, the initial velocity of the top plate 
along the z-axis was set as — le—7 m per step (minus represented that the velocity was 
along the negative direction of the z-axis). Since the velocity of the plate was fixed 
during the simulation process, the top plate would decline with a constant velocity. 
Additionally, no velocity was applied on the bottom plate. Consequently, the bottom 
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Fig. 3.2 Grid points monitored in the simulation process 


plate would be stable in the whole simulation process. With this method, the sample 
in the middle would be compressed by the top plate to simulate the UCS process. 
In the simulation process, the displacement along the z-axis of four grid points in 
the top plate was monitored (Fig. 3.2). 
Those four displacements were averaged as the uniaxial displacement of the 
sample, as shown in Eq. (3.1). 


s= Sal + Sa2 t Sa3 + Sa4 (3.1) 


where są: the axial displacement; and s,;: the axial displacement of the grid point i 
in the top plate. 

In UCS test, to calculate the Poisson’s ratio, the circumferential strain gauge was 
usually attached at the middle of the sample, to record the circumferential displace- 
ment of the sample. Then, in this simulation, the circumferential displacement of 
four grid points at the middle of the outside surface of the sample was monitored. 
Those four displacements were also averaged as the circumferential displacement of 
the sample, as shown in Eq. (3.2). 


= Sel + Sc2 t 53 + Sea (3.2) 


where se: the circumferential displacement of the sample; Sci: the circumferential 
displacement of the grid point i at the middle of outside surface of the sample. 
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Since the bottom plate was stabilised in the whole simulation process, the reaction 
force along the z-axis of the grid points at the bottom surface of bottom plate was 
summed. Then, this summed force was divided by the cross-section area of the 
sample, and the result was regarded as the uniaxial compressive stress of the sample. 

Then, the numerical UCS test was conducted. The stress-to-strain relationship of 
the sample is shown in Fig. 3.3. In the numerical simulation, the UCS of the sample 
was 55.87 MPa. 

In the experimental test, the UCS of the cement grout sample was 50.6 MPa, 
which was consistent with the numerical simulation results. 

It also shows that in numerical simulation, after the stress of the sample reached 
the peak, the circumferential strain decreased. To check the reason, the circumferen- 
tial displacement of a grid point was taken as an example and plotted. The grid point 
at the right middle of the sample (with an ID of 1430) was selected. The circumferen- 
tial displacement of the grid point with the calculating time step is shown in Fig. 3.4. 
The results show that before the calculating time step reached 6650, the circumfer- 
ential displacement of the grid point increased. This indicated that the diameter of 
the sample at the middle increased because of the compression. However, after the 
calculating time step reached 6650, the circumferential displacement of the grid point 
started decreasing. This indicated that after the sample reached the peak strength, 
the diameter of the sample in the middle also reached the peak. After the sample 
failed, the diameter of the sample in the middle started decreasing. This explained 
the phenomenon that after the sample reached the peak strength, the circumferential 
strain started decreasing. 

Based on the numerical simulation, the elastic modulus and the Poisson’s ratio can 
also be acquired. The results showed that the elastic modulus in numerical simulation 
was 13.47 GPa. The Poisson’s ratio can be calculated with Eq. (3.3) [31]: 


Ec 
p=- (3.3) 
Ea 
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where u: Poisson’s ratio; £.: the circumferential strain of the sample; and ¢,: the 
axial strain of the sample. 

Then, the Poisson’s ratio can be calculated as 0.26. Comparing the elastic modulus 
and Poisson’s ratio acquired from numerical simulation with those acquired from 
experimental tests, there was a close match between them, as shown in Table 3.2. This 
further confirmed that the input mechanical parameters were capable for simulating 
the cement grout. 


3.3.1.2 Numerical Direct Shear Test and Validation 


The input mechanical parameters used in the UCS simulation were then used in this 
section to simulate the direct shear process. Moosavi and Bawden [25] conducted 
experimental direct shear tests on cylindrical samples with a diameter of 100 mm and 
length of 200 mm. The test process was simulated in this section, and the simulation 
result was then validated with experimental results. 

In FLAC3D, a numerical cylindrical sample whose dimension was same as the 
experimental sample was established. The sample was composed of two sections, 
namely the top section and the bottom section. In the whole numerical model, there 
were 18,000 zones and 19,232 grid points. It should be mentioned that in the direct 
shear test, the shear failure plane is a pre-defined plane. Consequently, in numerical 
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i ee AN parameters Group Input parameter Value 
Interface Normal stiffness (GPa/m) 40 
Shear stiffness (GPa/m) 8.28 
Cohesion (MPa) 53 
Friction angle (°) 30 


simulation, an interface was created at the middle of the sample, namely between the 
top section and the bottom section, to simulate the shear failure plane. The interface 
was composed of 49,128 elements and 24,725 nodes. 

The constitutive model of the created interface was MC model, and the input 
parameters of the interface are tabulated in Table 3.3. 

The top section was fixed along the x-axis. For the bottom section, the external 
surface was fixed along the x-axis and y-axis. This was to simulate that in the exper- 
imental test, the bottom sample was restricted by the shear box and no movement 
could occur. As for the bottom surface of the sample, it was fully fixed to simulate 
that in the experiment, the whole sample was supported with a plate. A compressive 
normal stress of 2 MPa was applied on the top surface of the sample. Then, this 
numerical model was calculated until the unbalanced force ratio decreased to le—5. 
With this method, the initial normal stress equilibrium condition was acquired. The 
normal stress distribution in the sample is shown in Fig. 3.5. It is apparent to see 
that there was basically a uniform normal stress distribution in the sample. Under 
compression of the normal stress of 2 MPa, the normal stress in the sample was 
around 2 MPa. 

Before the direct shear test was performed, the displacement along the shearing 
direction of four grid points at the top surface of the sample was monitored. In the 
simulation process, the displacements of those four grid points were averaged as the 
shear displacement of the sample. 

It is more important to monitor the shear stress in the direct shear test process. To 
acquire this, in each step in the direct shear test process, each node on the interface 
was checked whether that node was still contact. If it was, the shear force of that 
node was recorded and summed. The accumulated shear force was the whole shear 
force of the interface and was then divided by the area of the interface. The result 
was regarded as the shear stress of the interface. 

Then, the direct shear test was performed. A shearing velocity of le—6 m per step 
was applied on the top section of the sample. This velocity was kept constant in the 
whole simulation process, to simulate the servo control. The direct shear test was 
conducted until the sample failed. 

The comparison between the numerical simulation result and the experimental 
test result is shown in Fig. 3.6. In numerical simulation, when the shear displace- 
ment increased to 0.94 mm, the shear stress of the interface increased to the peak 
of 6.48 MPa. After that, the shear stress decreased gradually. Also, there was a 
close match between numerical and experimental results. This confirmed that numer- 
ical simulation was capable in simulating the shear behaviour of the cement grout. 
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Fig. 3.5 Normal stress distribution in the sample 
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It should be mentioned that this comparison was compared with the final shear 
displacement being less than 2 mm. 

In the shearing process, the averaged normal stress of the interface was also 
monitored. The variation of the normal stress on the interface is shown in Fig. 3.7. 
Before the shear displacement increased to 0.94 mm, the normal stress of the interface 
increased slightly to the peak of 2.15 MPa. After that, the normal stress of the interface 
decreased gradually. However, in the whole process, the minimum normal stress of 
the interface was 1.99 MPa, which was only a little bit smaller than the maximum 
normal stress of the interface. This indicated that in the direct shear process, the 
normal stress of the interface was not constant although the applied normal stress 
was constant. This was probably due to the separation between the top sample and 
the bottom sample, and the frictional failure of the interface. 

After the direct shear test, a plot of the normal stress distribution on the interface 
is shown in Fig. 3.8. 
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The results showed that the normal stress at the left side of the interface was zero. 
This was due to the separation of the top and bottom samples. Consequently, the 
interface was non-contact around the left side, as shown in Fig. 3.9. 

Since it is more important to study the state of the interface that was contact, 
the stresses on the contacted interface were further studied. At the left side of the 
contacted interface, the normal stress was apparently much higher, with a maximum 
normal stress of 3.56 MPa. Moreover, the normal stress on the interface decreased 
towards the shearing direction. At the right side, the normal stress was much lower, 
with a minimum normal stress around 0.75 MPa. 
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Fig. 3.8 Normal stress distribution on the shear failure plane 
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Fig. 3.9 State of the interface in the direct shear test: a before the direct shear test; b after the direct 
shear test 


As for the shear stresses on the interface, its distribution is shown in Fig. 3.10. 
Around the left side of the whole interface, the shear stress was also zero, due to the 
non-contact of the interface. 

As for the contacted interface, the shear stress distribution was apparently not 
uniform. At the left side of the contacted interface, the shear stress was much higher, 
with a maximum of 7.36 MPa. Then, the shear stress decreased towards the shearing 
direction. At the right side, the shear stress was relatively lower, with a minimum 
around 6.0 MPa. This trend was basically consistent with the distribution of the 
normal stress on the interface. 
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Fig. 3.10 Shear stress distribution on the shear failure plane 


However, attention should be paid that Figs. 3.8 and 3.10 depicted the interface 
stress state after the test. In the middle shearing process, whether the distribution state 
was consistent with them was uncertain. And this will be studied in the following 
section. 

After the direct shear test, the plastic state of the sample is shown in Fig. 3.11. 
In this figure, “tension-n” means that the zone was experiencing tensile failure now. 
“tension-p” means that the zone previously experienced tensile failure. It is apparent 
to see that at the left side of the bottom sample, tensile failure occurred in the sample. 
Specifically, in the bottom sample, tensile failure occurred around the left side. In 
fact, this was consistent with the laboratory experimental test. Moosavi and Bawden 
[25] indicated that in the experimental test, tension crack occurred in the cement 
grout sample around the shear failure plane. 


3.3.2 Simulation of the Cement Grout with a Lower Strength 


Numerical simulation was also conducted on the same cement grout type while with 
a relative lower strength. First, the numerical UCS test was performed and compared 
with the experimental result, to confirm the proper input parameters. 


3.3.2.1 Numerical UCS Test and Validation 


The UCS test was simulated in FLAC3D. The simulation process was generally 
consistent with the procedures described above. Therefore, detailed simulation 
procedures were not described here. Only the simulation results were presented. 
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After the simulation, the stress-to-strain relationship of the sample was obtained. 
When the axial strain increased to 0.53%, the axial stress of the sample reached 
the peak of 46.09 MPa. In experimental test, the measured UCS was 40.3 MPa. 
Therefore, there was a close match between them. 

The elastic modulus and Poisson’s ratio were also obtained. In numerical simu- 
lation, the elastic modulus and Poisson’s ratio of the cement grout sample were 
8.33 GPa and 0.26. A comparison between experimental and numerical results 
is shown in Table 3.4. There was a good correlation between them. This further 
confirmed that the input parameters were robust for simulating the cement grout and 
can be used in the following direct shear test. 


3.3.2.2 Numerical Direct Shear Test and Validation 


Then, the numerical direct shear test was performed. And the input parameters of the 
sample were same as those used in the UCS simulation. A comparison between the 
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experimental and the numerical results is shown in Fig. 3.12. In numerical simulation, 
when the shear displacement increased to 1.41 mm, the shear stress of the interface 
reached the peak, which was 6.02 MPa. Following this, the shear stress reduced 
gradually with the shear displacement. The results also showed that there was a good 
correlation between experimental and numerical results. This further confirmed the 
effectiveness of numerical simulation in evaluating the shear behaviour of the cement 
grout. Attention was paid that in this comparison, the maximum shear displacement 
was less than 1.8 mm. 

Also, the variation of the normal stress on the shear failure plane was checked. In 
the simulation process, the normal stress on the shear failure plane was not constant. 
It had a slight increasing trend at the beginning and then decreased gradually. 


3.4 Analysis on the Shear Stress Distribution 


Figures 3.8 and 3.10 showed that after the direct shear test, the normal stress and 
shear stress distribution on the shear failure plane were not uniform. However, in the 
middle process of the direct shear test, the normal stress and shear stress distribution 
state were uncertain. Therefore, in this section, the distribution of the shear stress 
and normal stress on the shear failure plane was further studied. 

On the numerical shear stress versus shear displacement curve for the cement grout 
with a w/c ratio of 0.4, another three different points when the shear displacement 
increased to 0.6 mm, 0.94 mm and 1.2 mm, were selected. At those points, the shear 
stress and normal stress distribution on the shear failure plane was plotted. 

When the shear displacement was 0.6 mm, the shear stress and normal stress 
distribution on the shear failure plane was shown in Fig. 3.13. At this point, the 
shear stress distribution on the interface was generally uniform. On the contacted 
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interface, the maximum shear stress was 4.97 MPa while the minimum shear stress 
was 4.63 MPa. The relative difference between the maximum and the minimum shear 
stresses was only 0.34 MPa and the standard deviation was 0.09. Since this relative 
difference was too small, it can be neglected. 

As for the normal stress, on the contacted interface, there was a non-uniform 
distribution. Specifically, the maximum normal stress was around 2.96 MPa while 
the minimum normal stress was around 1.69 MPa. The relative difference between 
the maximum and the minimum normal stresses was 1.27 MPa, and the standard 
deviation was 0.29. 

When the shear displacement increased to 0.94 m, the distribution of the shear 
stress and normal stress is shown in Fig. 3.14. Currently, there was an apparent non- 
uniform distribution of the shear stress on the contacted interface. At the left side, the 
shear stress was relatively higher, with a maximum of 7.33 MPa. As for the right side, 
the minimum was around 6.20 MPa. The relative difference between the maximum 
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Fig. 3.13 Distribution of the shear stress and normal stress on the shear failure plane when the 
shear displacement arrived at 0.6 mm: a shear stress distribution; b normal stress distribution; c€ 
point selected on shear stress—shear displacement curve 
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shear stress and minimum shear stresses was 1.13 MPa, and the standard deviation 
was 0.28. 

The results also revealed that the non-uniform distribution of the normal stress 
became more apparent. On the contacted interface, the maximum normal stress was 
3.63 MPa, occurring at the left side while the minimum normal stress was 1.55 MPa, 
occurring at the right side. The relative difference between the maximum normal 
stress and minimum normal stress was 2.08 MPa, and the standard deviation was 
0.5. 

When the shear displacement reached 1.2 mm, the non-uniform shear stress distri- 
bution was still maintained, as shown in Fig. 3.15. On the contacted interface, the 
maximum shear stress was 7.39 MPa while the minimum shear stress was 6.18 MPa. 
The relative difference between the maximum shear stress and the minimum shear 
stress was 1.21 MPa, and the standard deviation was 0.29. 

As for the normal stress distribution, it was still non-uniform. At the left side 
of the contacted interface, the normal stress was much higher, with a maximum 


In Shear Stress 
7.3267E+00 
6.5940E+00 
5.8614E+00 
§.1287E+00 
4.3960E+00 
3.6634E+00 
2.9307E+00 
2.1980E+00 
1.4653E+00 
7.3267E-01 
0.0000E+00 


Norma! Stress 
3.6325E+00 
3.2692E+00 
2.9060E+00 
2.5427E+00 
2.1795E+00 
1.8162E+00 
1.4530E+00 
1.0897E+00 
7.2650E-01 
3.6325E-01 
0.0000E +00 


0.94 mm, 6.48 MPa 


Shear stress (MPa) 


—— Numerical results 
O Selected points 
T T 


0 1 f 
0 0.5 1 1.5 2 


Shear displacement (mm) 
(c) 


Fig. 3.14 Distribution of the shear stress and normal stress on the shear failure plane when the 
shear displacement arrived at 0.94 mm: a shear stress distribution; b normal stress distribution; c€ 
point selected on shear stress—shear displacement curve 
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Fig. 3.15 Distribution of the shear stress and normal stress on the shear failure plane when the 
shear displacement arrived at 1.2 mm: a shear stress distribution; b normal stress distribution; € 
point selected on shear stress—shear displacement curve 


of 3.61 MPa. At the right side, the normal stress was smaller, with a minimum of 
1.52 MPa. The relative difference between the maximum normal stress and minimum 
normal stress was 2.09 MPa, and the standard deviation was 0.5. 

The above three figures showed that at the beginning of the direct shear test, 
on the contacted area of the shear failure plane, the relative difference between the 
maximum shear stress and minimum shear stress was quite small. However, with the 
shear displacement increasing, the relative difference became larger. 

To show the variation of the relative difference of the shear stresses, in the 
calculating process, the relative difference between the maximum shear stress and 
the minimum shear stress was monitored, as shown in Fig. 3.16. When the shear 
displacement was smaller than 0.81 mm, the relative difference of shear stresses was 
quite small. This indicated that in this stage, the shear stress can be regarded that 
it uniformly distributed on the shear failure plane. However, after that, the relative 
difference of the shear stresses increased dramatically. When the shear displacement 
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Meanwhile, the relative difference between the maximum normal stress and 
minimum normal stress was monitored, as shown in Fig. 3.17. It shows that before the 
shear displacement increased to 0.97 mm, which was around the shear strength point, 
the relative difference increased monotonously. After that, the relative difference 


became also constant. 


point, the relative differ- 
ence of the shear stresses was maximal. After that, the relative difference became 
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3.5 Recommendation for Further Work 


Following the current study, the authors plan to continue the numerical study on the 
shear behaviour of cement grouts. Considering that the shear behaviour of cement 
grouts can also be studied with tri-axial tests, the authors may use the numerical 
method to evaluate the mechanical response of cement grouts under the tri-axial 
loading condition. Its purpose is to reveal the shear failure mechanism of cement 
grout during the tri-axial loading process. 


3.6 Conclusions 


Shear behaviour of the cement grout plays a significant role in determining the 
stability of systems and structures in civil and mining engineering. To better under- 
stand the shear behaviour of the cement grout, numerical direct shear tests were 
performed with the FLAC3D software. 

Cylindrical cement grout samples with two different strengths were established 
and numerically simulated. The numerical simulation results were compared with 
experimental results. There was a good correlation between them, confirming that 
numerical simulation can be used to evaluate the shear behaviour of the cement grout. 

The numerical direct shear test results showed that in the shearing process, along 
the shear failure plane, the shear stress at each point was not identical. There was a 
shear stress distribution on the contacted shear failure plane. At the beginning of the 
direct shear test, the relative difference of the shear stresses was quite small. Then, 
in this stage, it can be assumed that there was a uniform shear stress distribution on 
the contacted shear failure plane. However, once the shear displacement increased to 
around the shear strength point, the relative difference of the shear stresses became 
much larger. In this stage, the shear stress was non-uniformly distributed on the 
contacted shear failure plane. After that, the relative difference of the shear stress 
became stable. 

Additionally, although the applied normal stress on the sample was constant, 
the normal stress on the contacted shear failure plane was variable. Before the shear 
strength point, the normal stress increased slightly. After that, the normal stress grad- 
ually decreased. There was also a distribution of the normal stress on the contacted 
shear failure plane. This normal stress distribution was not apparent at the begin- 
ning of the test. However, when the shear displacement increased to around the 
shear strength point, it became much obvious. Moreover, after that, this non-uniform 
normal stress distribution became stable. 
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Chapter 4 A) 
Analytical Modelling to Study the Load get 
Transfer Performance of Fully Grouted 

Rock Bolts 


4.1 Introduction 


A rock bolt is a rock tendon which can be either a rod or a strand, installed in a 
borehole to reinforce the surround rock mass and soil [1]. Rock bolts have been used 
in civil and mining for a long time [2]. Numerous filed applications and computer 
simulation have proved that rock bolts play a significant role in improving the internal 
strength of rock mass and enhancing the safety of underground excavations [3—6]. 

The rock bolt reinforcement system is basically composed of four components: the 
surrounding confining medium, e.g. rocks or soils, the rock reinforcement element, 
the internal fixture and the external fixture [7—9]. It is found that all these four 
components have an influence on the performance of a reinforcement system. 

To better understand the performance of the reinforcement system, Thompson 
et al. [10] classified the rock reinforcement system as three types: Continuously 
Mechanically Coupled (CMC), Continuously Frictionally Coupled and Discrete 
Mechanically or Frictionally Coupled. This study deals with fully grouted rock bolts, 
and based on this classification, fully grouted rock bolts belong to the CMC system. 

For the CMC system, cementitious grout or polyester resin can be used to bond 
rock bolts with surrounding confining medium [11-13]. Via the bonding material, 
the load from exterior unstable rock mass can be transferred to interior stable rock 
mass. It should be mentioned that in the in situ applications, rock bolts can encounter 
axial, shear, combined axial and shear loading conditions [14]. However, this paper 
focused on the axial performance of rock bolts. 

Failure of the rock bolting system has different modes. Jeremic and Delaire [15] 
summarised that failure of the rock bolting system can be classified as bond failure of 
the bolt/grout interface, bond failure of the grout/rock interface, failure of the grout 
annulus and failure in the surrounding rock mass. Then, Hutchinson and Diederichs 
[16] added that rupture of the rock tendon also occurred in the in situ scenario. 
Nevertheless, among these failure modes, bond failure of the bolt/grout interface and 
the grout/rock interface are more common to encounter [17]. 
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Jalalifar [18] stated that when fully grouted rock bolts were loaded axially, the 
shear stress along the bolt/grout interface was larger than the shear stress along 
the grout/rock interface because the bolt/grout interface had a smaller contact area. 
When the grout and the rock had similar strength, failure was likely to occur along 
the bolt/grout interface. Previous experimental pull-out tests also showed that bond 
failure of the bolt/grout interface occurred in the tests [19, 20]. Therefore, this paper 
deals with the bond failure of the bolt/grout interface. 

Gambarova [21] stated that the bonding performance of the bolt/grout interface 
mainly relies on chemical adhesion, mechanical interlock and friction. Chemical 
adhesion plays a negligible role on the bonding performance and can be easily 
removed after a small relative movement between the bolt and grout [22]. 

To understand the load transfer performance of rock bolts, extensive research has 
been conducted via experimental pull-out tests [23-25]. Compared with that, much 
less work has been conducted with analytical modelling. Hawkes and Evans [26] 
conducted the pioneering work in analytically modelling the load transfer of rock 
bolts. Specifically, an analytical model was proposed to evaluate the load distribution 
along a fully grouted rock bolt. However, this assumption was only valid in the 
coupling stage. In the modelling process, a linear relationship between the axial stress 
in the bolt and the shear stress along the bolt/grout interface was assumed. However, 
this assumption was credible when the bolt was fully coupled with the surrounding 
medium. With the axial load in the bolt increasing, debonding occurred along the 
bolt/grout interface and this assumption may not be valid again. Farmer [27] deduced 
an analytical model to study the load distribution along a rock bolt. In this analytical 
model, it was assumed that the shear stress along the bolt/grout interface was linear 
with the grout shear modulus. Therefore, this model was only applicable when the 
bolt/grout interface underwent elastic behaviour. Once debonding occurred along the 
bolt/grout interface, the credibility of this analytical model decreased. Aydan et al. 
[28] deduced an analytical solution to study the load transfer of rock bolts. However, 
no experimental work was conducted to validate the credibility of the model. Rajaie 
[29] studied the load transfer of cable bolts and an energy approach was adopted to 
analyse the loading process of cable bolts. There was a good correlation between 
laboratory results and analytical results. However, this model mainly focused on the 
peak load calculation and no equation was deduced to calculate the load variation 
with displacement in the pull-out process. Therefore, this model cannot be used to 
study the load—displacement relationship of cable bolts. 

Yazici and Kaiser [30] proposed a bond strength model to study the load transfer 
of plain cable bolts. Later, this model was modified to study the influence of stress 
change on load transfer of cable bolts [31]. This model treated the bolt/grout inter- 
face as a zigzag geometry, similar to a rock joint. For the rock joint, strength of 
the rock is usually weak along two sides. However, for the cable bolting system, 
strength of the steel is much higher along the bolt side. Following that, Hyett et al. 
[32] developed a constitutive model to study the load—displacement performance of 
cable bolts. Two different boundary conditions, constant normal stress and constant 
normal stiffness, were considered. This model assumed that the shear stress along 
the bolt/grout interface was uniform. However, experimental work proved that when 
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the bolt had a long embedment length, the shear stress along the bolt/grout interface 
was not uniform [22]. Therefore, this model was applicable to calculate the load— 
displacement relationship for short-embedded bolts. However, it may not be suitable 
for long-embedded bolts. A load distribution model was proposed by Hyett et al. 
[33]. This model can study the load distribution in the embedment section and the 
anchor section. Moosavi [34] developed an analytical model for modified cable bolts. 
However, when modified cable bolts were used to reinforce rock masses, the modi- 
fied geometry can be installed at different positions in the embedment section. And, 
the influence of the modified geometry position on the load transfer performance of 
cable bolts was not considered. Li and Stillborg [35] developed an analytical model 
to study the load distribution of rock bolts. The bolt/grout interface was assumed to 
be composed of four sections: completely debonded, partially debonded, softening 
and elastic. 

Cai et al. [36] used a modified shear-lag model to study the load transfer of rock 
bolts. However, the softening behaviour of the bolt/grout interface was neglected. 
Xiao and Chen [37] assumed that the bolt/grout interface underwent elastic, softening 
and debonding stages and developed a model for rock bolts. Ren et al. [38] deduced 
the load transfer process of rock bolts with an analytical model. Experimental work 
was used to validate the credibility of this analytical model. However, when deducing 
the analytical model, it was assumed that the shear slippage of the bolt/grout inter- 
face was equal to the axial displacement of the rock bolt. Then, the displacement of 
the surrounding confining medium was neglected. Consequently, this model cannot 
be used to study the influence of the surrounding confining medium on the perfor- 
mance of rock bolts. Blanco Martin et al. [39] proposed an analytical solution to 
evaluate the load—displacement performance of rock bolts. Jalalifar [40] evaluated 
the performance of rock bolts installed in an elastoplastic rock mass. In this model, 
it was assumed that the shear force in the grout annulus had a linear relationship 
with the shear slippage. Therefore, this model was only valid when the bolt/grout 
interface underwent elastic deformation. After the bolt/grout interface started plastic 
deformation, this model was not valid again. Cao et al. [41] studied the influence of 
rib geometry on the bonding capacity of rock bolts. It was predicted that rock bolts 
with small rib face angles were more effective in hard rock conditions. 

Ma et al. [42] used a nonlinear equation to depict the bond-slip behaviour of the 
bolt/grout interface and studied the load transfer of rock bolts. Later, the model was 
modified to study the influence of free-end slip on the performance of rock bolts 
[43]. However, in this model, an empirical equation was used to depict the bond- 
slip relationship of the bolt/grout interface. Therefore, the properties of the rock 
bolting system may modify the coefficients in the bond-slip model. Consequently, 
this model cannot be used to predict the influence of rock mass confinement on 
the performance of rock bolts. Chen et al. [13] used an analytical model to study 
the axial performance of cable bolts. However, the model was only appropriate for 
laboratory pull-out scenario. Li et al. [44] developed an analytical model to study 
the load transfer of modified cable bolts. Cable bolts with short embedment length 
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were used to validate the model. However, in situ tests proved that the installed cable 
bolts can be very long for the field scenario [45]. And whether this model can be 
used to study the performance of cable bolts with long embedment length was not 
discussed. 

The literature review shows that for the rock bolting system, when the shear 
behaviour of the bolt/grout interface was depicted with a bond-slip equation, different 
analytical models can be used, such as the exponential model [42] or the tri-linear 
model [46]. Among them, the tri-linear model considers the elastic, softening and 
debonding behaviour of the bolt/grout interface and thus is more commonly used. 
However, when the tri-linear model was incorporated into the rock bolting system 
to deduce the load—displacement relationship of rock bolts, the involvement of the 
surrounding confining medium was usually not considered. 

It should be mentioned that Vaculik et al. [47] conducted an analytical study to 
evaluate the load-slip relationship of the FRP-to-substrate joints with the bond-slip 
model. A framework was proposed to depict the performance of FRP-to-substrate 
joints. However, for the rock bolting system, cylindrical rock bolts are always used 
to reinforce the excavated rock masses, such as tunnels, roadways and stopes [48, 
49]. Therefore, this paper applies the method proposed by Vaculik et al. [47] into the 
fully grouted rock bolting system and further deduced the specific analytical solution 
for fully grouted rock bolts. Additionally, they evaluated the performance of FRP- 
to-substrate joints with the embedment length ranging from short to long. However, 
for fully grouted rock bolts, the embedment length is usually extremely long [45, 
50]. For example, Cao et al. [51] summarised that the embedment length for a rock 
bolt can be around 3 m. Furthermore, when the cable bolt is used, the embedment 
length may range from 3 to 15 m. Considering this, this paper especially focused on 
the long embedment length scenario for rock bolts. 

First the relationship between the shear stress of the bolt/grout interface and the 
shear slippage was simulated with a tri-linear model. Then, the pull-out process of 
rock bolts was analysed with five stages: elastic, elastic-softening, elastic-softening- 
debonding, softening-debonding and debonding stages. After that, the analytical 
model was validated with experimental pull-out tests. Last, a parametric study was 
conducted to evaluate the influence of elastic modulus of the confining medium, 
the shear strength of the bolt/grout interface and the residual shear strength of the 
bolt/grout interface on the load transfer performance of rock bolts. 


4.2 Analytical Modelling 


4.2.1 Bond-Slip Model 


A tri-linear model was used to depict the bond-slip behaviour of the bolt/grout inter- 
face. The tri-linear model is composed of three parts: ascending section, descending 
section and horizontal section [13, 38], as shown in Fig. 4.1. Specifically, with the 
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Fig. 4.1 A tri-linear model T 

[13, 38, 47] 
Tp ——— 
Tr 


Ôr ô 


shear slippage of the bolt/grout interface increasing from 0 to ôp, the shear stress 
increases linearly from zero to the shear strength of the bolt/grout interface. Then, 
with ò increasing to §,, the shear stress decreases linearly to the residual shear strength 
of the bolt/grout interface. After 5,, the shear stress keeps constant and equals t+. 
Therefore, the tri-linear model can be expressed as Eq. (4.la—4.1c) [13, 38]: 


c= 28 (0<8 <4) (4.1a) 


r= P —*§ (5) <5 <6) (4.1b) 


T=T, (6 > 6) (4.1c) 


where Tp: shear strength of the bolt/grout interface; t,: residual shear strength of the 
bolt/grout interface; and ô: shear slippage when the residual shear strength of the 
bolt/grout interface reaches. 


4.2.2 Governing Equation 


In previous research, the shear slippage of the bolt/grout interface was usually simpli- 
fied as the axial displacement of the rock bolt [38, 42]. However, previous experi- 
mental work proved that the shear slippage of the bolt/grout interface was different 
from the axial displacement of the rock bolt because of the displacement of the 
confining medium [52]. Furthermore, Hyett et al. [33] found that the displacement 
of the confining medium had an apparent influence on the load transfer performance 
of rock bolts. Therefore, neglecting the displacement of the confining medium in 
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the analytical model could not properly reflect the load transfer performance of rock 
bolts. Considering this issue, in this paper, the displacement of the confining medium 
was included. 

Since this study aims at studying the bond failure of the bolt/grout interface, the 
grout/medium interface was assumed to be intact in the pull-out process of rock 
bolts. Then, the shear slippage of the bolt/grout interface was defined as the relative 
movement between the rock bolt and the confining medium [53, 54]: 


d(x) = Ub (x) — Um(x) (4.2) 


where 5(x): shear slippage of the bolt/grout interface at the position of x; u(x): rock 
bolt displacement at the position of x; and u(x): confining medium displacement at 
the position of x. 

Differentiating Eq. (4.2) twice with respect to x leads to: 


Sa) dux) dunl) 
agl 2 (4.3) 
x dx dx 


The axial stress in the rock bolt can be expressed with the displacement of the 
rock bolt. As for the confining medium, this is also applicable. Then, the following 
equations can be acquired [40]: 


olx) = Es mo) (4.4) 
X 
o duma) 
Om(xX) = Em (4.5) 
dx 


where op(x): axial stress in the rock bolt at the position of x; Em: elastic modulus of 
the confining medium; and o,,(x): axial stress applied on the confining medium. 
Substituting Eqs. (4.4) and (4.5) into (4.3) leads to: 


PS) 1 dax) 1 dom(x) 
dx? Ep dx Em dx 


(4.6) 


According to the mechanical equilibrium, the relationship between the axial stress 
in the bolt and the axial stress applied on the confining medium can be expressed as 
[B3]: 


2 
Tri 


ee E 4.7 
o; i (4.7) 


where Apm: cross-section area of the confining medium; and rẹ: rock bolt radius. 
Substituting Eq. (4.7) into Eq. (4.6) leads to: 
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Fig. 4.2 The cross-section confining medium 
view of a rock bolting system a 


o, <—— bolt — > o;+do, 


BE- grout column 


sax) (1 mre \ do (x) (48) 
dx?  \Ep EmAm/) dx 


The relationship between the axial stress in the rock bolt and the shear stress along 
the bolt/grout interface can be seen in Fig. 4.2. 
Then, the following equation can be acquired [42]: 


rp dop (x) 
= 4.9 
t(x) =e (4.9) 
where t(x): shear stress of the bolt/grout interface at the position of x. 
Substituting Eq. (4.9) into Eq. (4.8) leads to: 
s(x) 2/1 mre (x) =0 (4.10) 
T(x) = . 
dx? Fb Ey EmAm 
It is defined that: 
2/1 2 
ees E (4.11) 
rp Ep EmAm 
Then, Eq. (4.10) can be written as: 
dS) a 
ro oe à t(x) =0 (4.12) 


Equation (4.12) is the governing equation of this analytical model, and it can be 
solved with the pre-defined shear stress equations. 


4.2.3 Pull-Out Stages of Rock Bolts 


According to previous research, the pull-out stages of the bolt/grout interface can 
be classified as elastic stage, elastic-softening stage, elastic-softening-debonding 
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stage, softening-debonding stage and debonding stage [38, 53, 56]. Therefore, in the 
following contents, the pull-out process of a rock bolt was analysed with those five 
stages. 

4.2.3.1 Elastic Stage 

When the pull-out load is small, the bolt/grout interface undergoes elastic deforma- 
tion. And the full bolt/grout interface is elastic, as shown in Fig. 4.3. 


Then, Eq. (4.1a) can be substituted into Eq. (4.12), leading to: 


d? (x) D 


P42 
d(x) =0 4.13 
ae 7 P80) (4.13) 
It is defined that: 
= fi (4.14) 
p 


Then, Eq. (4.13) can be written as: 


d8 (x) 


Ted xax) =0 (4.15) 


Equation (4.15) can be solved with boundary conditions. First, the axial stress in 
the bolt should be deduced. According to Eq. (4.2), the following equation can be 
acquired: 


dd(x) _ dup(x) dum (x) 
dx — dx dx 


(4.16) 
It can be further written as: 


Confining medium 


Free end a Grout Loaded end: 
rs 


Db 
ro] 
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Fig. 4.3 The full bolt/grout interface is elastic in the elastic stage 
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d3(x) w) Om) 


4.17 
dx Ey Em ( ) 
Substituting Eq. (4.7) into Eq. (4.17) leads to: 
dd(x) 1 4: mre (x) (4.18) 
= Opl X m 
dx Ep EmAm b 


Substituting Eq. (4.11) into Eq. (4.18), the axial stress in the bolt can be written 
as: 


(x) 2 déd(x) (4.19) 
w(x) = —> ; 
b pÀ? dx 
The boundary conditions are: 
olx =0)=0 (4.20) 
F 
oa = L) = — (4.21) 
nre 
Combining Eqs. (4.19), (4.20) and (4.21) with Eq. (4.15) leads to: 
Fàiôp 
5(x) = ————_———- cosh (àx) (4.22) 
27r Tp Sinh(A; L) 
= ————_ sinh(a 4.23 
= SE) ay 
(x) ee h(x) (4.24) 
T(x) = ——————__ cos x ; 
mrp sinh (A1 L) : 


It is defined that at the loaded end, the shear slippage of the bolt/grout interface 
is equal to the pull-out displacement of the rock bolt [57]. Substituting x = L into 
Eq. (4.22), the pull-out load of the rock bolt can be expressed as: 


2 tanh(A,L 
ee) (4.25) 
À1ôp 
At the end of this stage, the shear stress at the position of x = L increases to the 
shear strength of the bolt/grout interface. Therefore, substituting t(x = L) = Tp into 
Eq. (4.24) leads to: 


27 rpTp tanh(A; L) 


Foop = 4 


(4.26) 
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where F'sof: pull-out load when the bolt/grout interface starts softening. 


4.2.3.2 Elastic-Softening Stage 


In the elastic-softening stage, the full bolt/grout interface is composed of elastic 
section and softening section, as shown in Fig. 4.4. 

It is defined that a, is the softening length. In the elastic-softening stage, the 
interface in the domain [0, L — ag] is still elastic. A boundary condition is Eq. (4.20) 
and the other boundary condition is: 


r(x = L—a,) =t (4.27) 


Combining Eqs. (4.9), (4.19), (4.20) and (4.27) with Eq. (4.15) leads to: 


— P 
80) = aa y oma) (4.28) 
As T sinh(A1x) (4.29) 
rpA, cosh(A;(L — as)) 
Z Tp 
t(x) = SCE a) cosh(A;x) (4.30) 


Substituting x = L — as into Eq. (4.29), the axial stress can be acquired: 
2Tp 
(x = L — as) = —— tanh (à; (L — as)) (4.31) 
rpAt 


In the domain [L — a,, L], the bolt/grout interface softens. Substituting Eq. (4.1b) 
into Eq. (4.12) leads to: 
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Fig. 4.4 The full bolt/grout interface is composed of elastic section and softening section in the 
elastic-softening stage 
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PE) 5% T 3 Tpdr — Trôp 
À 6 À =0 4.32 
rs F a ) S, — ô» Sa 
It is defined that: 
Tp — T, 
Ma 2 7}? 4. 

e 3, (4.33) 

Therefore, Eq. (4.32) can be further written as: 

PS) a 2 Tpôr — qô 

158x) — 13 2—2 =0 4.34 
dx? PAPE Tp — Tr man 


The boundary conditions are Eqs. (4.27) and (4.31). Combining Eqs. (4.9), (4.19), 
(4.27) and (4.31) with Eq. (4.34) leads to: 


d(x) = Aito ~ tanh(A,(L — a;)) sin(A2(x + as — L)) 


2 Tpôr — Trdp 


es 
= =~; cos(à2(L — as — x)) + (4.35) 
Aa Tp — Tr 
2t (1 
(x) = (= tanh(A;(L — a,)) cos(A2(L — as — x)) 
++ sin(à2(x + as — D») (4.36) 
2 


t(x) = o( 2 tanh(A;(L — a,)) sin(A2(L — as — x)) + cos(à2 (L — a, — ») 
(4.37) 


Substituting x = L and o (x = L) = , into Eqs. (4.35) and (4.36), respectively, 
the pull-out displacement and pull-out load of the rock bolt can be acquired: 


1 1 
F= 2am (= tanh(A,(L — ag)) cos(Azas) + z sin(aa.)) (4.38) 
1 


126, Tpdr — TrSp 


A416 
Up = = us tanh(A; (L — ag)) sin(Azas) — 42 —1  cos(A2ds) + (4.39) 
2 


Pp r 
At the end of this stage, the shear stress at the position of x = L reduces to T,. 


Therefore, substituting t(x = L) = T, into Eq. (4.37) leads to: 


À 7 
cos(A2a,) — = tanh(A,(L — a;)) sin(a2a;) = = (4.40) 
1 P 
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For an infinite embedment length, tanh(A (L — a,)) is equal to 1. Then, the softening 
length when the elastic-softening stage ends can be acquired: 


Fi TAY 
arcsin (4.41) 


Ài 
Vat tay Try At +13 


where ası: the softening length when the elastic-softening stage ends. 


arcsin 


as) = = 


1 
A2 


4.2.3.3 Elastic-Softening-Debonding Stage 


After t(x = L) reduces to T,, the bolt/grout interface enters the elastic-softening- 
debonding stage. In this stage, the full bolt/grout interface is composed of elastic 
section, softening section and debonding section, as shown in Fig. 4.5. 

Assuming that the debonding length is ag, the bolt/grout interface in the domain 
[0, L — ag — as] is still elastic. Therefore, Eqs. (4.28), (4.29) and (4.30) are still valid 
if L is replaced with L — ag: 


— 51 
8) = Sa aay comme (4.42) 
ete. 2p sinh(A1x) (4.43) 
rpA, cosh(A,(L — dag — as)) 
Tp 


t(x) = 


cosh(A;(L — ag — as)) cosh(A1x) (4.44) 


As for the domain [L — ag — as, L — aq], the bolt/grout interface softens. Therefore, 
Eqs. (4.35), (4.36) and (4.37) are still valid if L is replaced with (L — aa): 
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Fig. 4.5 The full bolt/grout interface is composed of elastic section, softening section and 
debonding section in the elastic-softening-debonding stage 
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A18 
8(x) = a tanh(A;(L — aq — as)) sin(à2 (x + ag + a; — L)) 
2 


does Tpôr — Trdp 
— ~ cos(Aa(L — ag— as — x)) + (4.45) 
25 Tp — Tr 
2t, (1 
(x) = — a tanh(àı (L — dg — ds)) cos(à2 (L — dg — as — x)) 
Fb 1 
1 
E sin(à2(x + ag + a; — 1») (4.46) 
2 
ha : 
t(x) = (2 tanh(A;(L — aq — as)) sin (à2 (L — dg — ds — x)) 
1 
+cos(à2(L — ag — as — x))) (4.47) 
Substituting t(L — aa) = T, into Eq. (4.47) leads to: 
do r Tr 
cos(àz2as) — 1. tanh(àı (L — dag — ds)) sin(A2as) — = =0 (4.48) 
1 p 


Equation (4.48) depicts the relationship between ag and a,. It shows that in this 
stage, the softening length keeps varying with the debonding length increasing. 
Combining Eqs. (4.45) and (4.48) together, the shear slippage of the bolt/grout 
interface at the position of x = L — ag can be acquired: 


ôx = L — aa) = ô, (4.49) 


Substituting x = L — ag into Eq. (4.46), the axial stress can be acquired: 
2t (1 
(x = L — aa) = — | — tanh (à; (L — ag — ag)) cos (à2as) 
ro Ai 
l. 
+— sin(A2a;) (4.50) 
A2 


For the domain [L — ag, L], the bolt/grout interface debonds and the shear stress 
equals T;: 


TX) =T; (4.51) 


Also, Eq. (4.1c) can be substituted to Eq. (4.12), leading to: 


d8 (x) -x 


o T (4.52) 


Substituting Eqs. (4.19), (4.49) and (4.50) into Eq. (4.52) leads to: 
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l2 2 2 TA? 
d(x) = ôr + =" t(x — (L — aa) ) + Ea tanh(A;(L — aq — ds)) cos(A2as) 
1 


tae Z 2 
E sin(Aza,) — à T, (L — aa) | (x + aa — L) (4.53) 
2 


oe / 1 1 
OE 2h tatty (L — aq — a) cosa + — sin(aa)) 
ro \Aq d2 
2T: 
+ æ+- I) (4.54) 
b 


Substituting x = L and op(x = L) = F into Eqs. (4.53) and (4.54), the pull-out 


displacement and load of the rock bolt can be acquired: 


ai Te 


Up = ôr + 2 


TA Tpà? . 

+ ca tanh(A;(L — dag — ds)) cos(Aza,) + Sa sin(A2as) Jaa 
1 2 

(4.55) 


F= 2an (2 tanh(à1 (L — aa — ad) cos(àzas) + = sin(Aza,) + na) (4.56) 
1 2 


The maximum pull-out load occurs in this stage and can be acquired by solving 
Eqs. (4.48) and (4.56) together. 

At the end of this stage, the elastic length disappears. Therefore, L — ag = as. 
Substituting it into Eq. (4.48) leads to: 


1 Te 
ag = T arccosį| — (4.57) 


where as2: softening length when the elastic-softening-debonding stage ends. 


4.2.3.4 Softening-Debonding Stage 


After the elastic length disappears, the bolt/grout interface enters softening- 
debonding stage. In the domain [0, L — ag], the bolt/grout interface softens. Substi- 
tuting Eqs. (4.9), (4.19) and (4.20) and t(x = L — aa) = T, into Eq. (4.34) leads 
to: 


T(d5 — ôr Oy — Td 

( P ) cos(à2x) + Tpôr — trp 

(tp — Tr) cos(à2 (L — aa)) Tp — Tr 
Qt 


op(x) = can. oa) sin(A2x) (4.59) 


8(x) = (4.58) 
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T(x) = cos(A2x) (4.60) 


cos(A2(L — aq)) 
Substituting x = L — ag into Eqs. (4.58) and (4.59) leads to: 


ô(x = L — aa) = ô (4.61) 


PA 
(x = L — aa) = ER tan(A2(L — aa)) (4.62) 
pA2 


For the domain [L — ag, L], the bolt/grout interface debonds. Substituting Eqs. (4.19), 
(4.61) and (4.62) into Eq. (4.52) leads to: 


AT 


5 (x? —-(L- aa)’) 


5G) = ô; + 


2 
+ (= tan(A2(L — aa)) — X? (L — aa) (x +aa— L) (4.63) 
wh 


2t, ( tan(A2(L — aa)) 
(x) = ( 


+x+da- L) (4.64) 
A2 


Substituting x = L and o4(x = L) = F- into Eqs. (4.63) and (4.64), the pull-out 


5 
wr, 


displacement and load of the rock bolt can be acquired: 


tan(A.(L — 
ii tana S y Ua 2) (4.65) 
2 do 
tan(A.(L — 
F= dary, ( AOE =) + aa) (4.66) 
2 


At the end of this stage, the debonding length equals the full embedment length. 
Therefore, substituting ag = L into Eqs. (4.65) and (4.66), the pull-out displacement 
and load of the rock bolt at the end of the softening-debonding stage can be acquired: 


XGL? 
Ubo = ôr + 2 (4.67) 
Fo = 2TFpT;L (4.68) 


where upo: pull-out displacement of the rock bolt at the end of the softening- 
debonding stage; and Fọ: pull-out load of the rock bolt at the end of the 
softening-debonding stage. 
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4.2.3.5 Debonding Stage 

After the softening length disappears, the bolt/grout interface enters debonding stage. 
Since the full bolt/grout interface undergoes debonding, the relationship between the 
pull-out load and displacement can be expressed as: 


Fo — F = 27 rpt, (Up — uno) (4.69) 


Substituting Eqs. (4.67) and (4.68) into Eq. (4.69) leads to: 


XGL? 
F =2nrpt,| L+ ôr + z T (4.70) 


4.3 Calibration of the Input Parameters 


In this analytical model, the mechanical properties of the rock bolt and the confining 
medium, such as rp, E» and Em, can be acquired based on the rock bolt pull-out 
scenario. Specifically, r, and E, can be acquired from the rock bolt specification 
provided by the bolt manufacturer. Em can be acquired by conducting uniaxial 
compressive strength tests on samples cored from the confining medium. 

As for the mechanical properties of the bond-slip model, they can be calibrated 
from the experimental load—displacement results. Figure 4.6 shows a typical load— 
displacement curve of rock bolts. 

In the initial stage of the load—displacement curve, the pull-out load of the rock bolt 
generally had a linear relationship with the pull-out displacement. Therefore, Point 
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1 represents the point where the linear relationship of the pull-out load and displace- 
ment ends in the load-displacement curve. Then, 6, equals the pull-out displacement 
at Point 1. Also, F’so¢ equals the pull-out load at Point 1. Substituting F'sof and 6, into 
Eq. (4.26), Tp can be acquired. 

When the rock bolt reaches the maximum pull-out load, that point can be selected 
as Point 2, as shown in Fig. 4.6. Then, ô; can be selected in the range between the 
pull-out displacement at Point | and the pull-out displacement at Point 2. Then, t, can 
be calibrated until the maximum analytical pull-out load fits well with the maximum 
experimental pull-out load. 


4.4 Validation of the Analytical Model 


In this section, two experimental pull-out tests were used to validate the credibility 
of this analytical model. 


4.4.1 Validation with a Pull-Out Test 


Baiet al. [58] conducted experimental in situ pull-out tests on rock bolts. Specifically, 
rock bolts with a radius of 14 mm were tested and the embedment length was 3 m. 
The surrounding rock mass had an elastic modulus of 6 GPa. 

The mechanical properties of the rock bolt and the confining medium used in this 
analytical model are tabulated in Table 4.1. 

As for the input parameters for the bond-slip model, they can be calibrated from 
experimental pull-out results, following the methods indicated in Sect. 4.3. Then, the 
parameters tabulated in Table 4.2 were used in this paper to analytically model the 
load transfer performance of rock bolts. 

A comparison between the experimental result and the analytical result is shown 
in Fig. 4.7. There is a close match between the experimental result and the analytical 
result, confirming the credibility of this analytical model. 


Table 4.1 Mechanical properties of the bolt and the confining medium 
rp (mm) Ep (GPa) L (m) Em (GPa) Am (m?) 
14 51 3 6 1 


Table 4.2 Input parameters 


MP: ô MP ô 
for the bond-slip model Tp (MPa) pam Te (MPa) r (mm) 


2.2 3.57 1 8.91 
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4.4.2 Validation with the Other Pull-Out Test 


Stillborg [59] conducted a number of pull-out tests in the mine site to study the load 
transfer performance of cable bolts. Cable bolts with a long embedment length of 10 m 
were tested. The cable bolts had a radius of 19 mm and elastic modulus of 83 GPa. The 
surrounding rock mass had an elastic modulus of 57 GPa. The mechanical properties 
of the rock bolt and the confining medium used in this analytical model are tabulated 
in Table 4.3. 

As for the input parameters for the bond-slip model, they can be calibrated from 
the experimental pull-out results. In this paper, the parameters tabulated in Table 4.4 
were used to calculate the load—displacement relationship of the cable bolt. 

The comparison between the experimental result and the analytical result is shown 
in Fig. 4.8. It shows that the analytical result agrees well with the experimental result, 
further confirming that this analytical model can be used to simulate the load transfer 
performance of rock bolts. 


Table 4.3 Mechanical properties of the cable bolt and the confining medium 
rp (mm) Ey (GPa) L (m) Em (GPa) Am (m?) 
19 83 10 57 1 


Table 4.4 Input parameters 


for the bond-slip model i oy (rami) Te MPa) ô: (mmi) 


1.34 10.37 0.47 35.02 
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4.5 Parametric Study 


A parametric study was conducted to evaluate the input parameters on the load 
transfer performance of rock bolts. Specifically, the influence of the elastic modulus 
of the confining medium, the shear strength of the bolt/grout interface and the shear 
slippage where shear strength of the bolt/grout interface reaches on the load transfer 
performance of rock bolts were studied. 


4.5.1 Elastic Modulus of the Confining Medium 


To study the influence of the elastic modulus of the confining medium on the load 
transfer performance of rock bolts, an analytical rock bolt pull-out scenario was 
created. A rock bolt with a radius of 10 mm was tested. The elastic modulus of the 
rock bolt and the embedment length were 196 GPa and 1.5 m. Then, the mechanical 
properties of the bolt and the confining medium are tabulated in Table 4.5. As for 
the input parameters for the bond-slip model, they are tabulated in Table 4.6. 

The elastic modulus of the confining medium was changed from 10 to 90 MPa 
and the load—displacement performance of rock bolts was calculated, as shown in 
Fig. 4.9. It shows that the elastic modulus of the confining medium had an apparent 


Table 4.5 Mechanical properties of the bolt and the confining medium 


rp (mm) E, (GPa) L(m) Em (MPa) Am (m?) 
10 196 1.5 10 1 

50 

90 
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Table 4.6 Input parameters 


MP. ô MP ô 
for the bond-slip model i) Spn a 


2 1.5 0.8 39 


influence on the load transfer performance of rock bolts. The initial stiffness and 
the peak load of the rock bolting system increased with the elastic modulus of the 
confining medium. 

This finding agrees well with previous research. Hyett et al. [60] conducted labora- 
tory and in situ pull-out tests on plain cable bolts. In the laboratory, PVC, aluminium 
and steel pipes were used to confine cable bolts. The results showed that the peak 
load of cable bolts increased with the elastic modulus of the confining medium. As 
for the in situ case, cable bolts were pulled from shale, limestone and granite. The 
results showed that when the confining medium changed from shale to granite, the 
peak load of cable bolts increased apparently. 

Figure 4.9 also shows that with the elastic modulus of the confining medium 
increasing, the increasing rate of the peak load was not constant. To further study the 
influence of the elastic modulus of the confining medium on the peak load, a series 
of analytical pull-out tests was conducted. The relationship between the rock bolt 
peak load and the elastic modulus of the confining medium was acquired, as shown 
in Fig. 4.10. 

It shows that for this scenario, before the elastic modulus of the confining medium 
increased to 280 MPa, the peak load increases rapidly. After that, the peak load is still 
increasing but with a much lower slope. This indicates that when the elastic modulus 
of the confining medium was higher than a critical limit, it had marginal effect on 
the peak load of rock bolts. 


Fig. 4.9 Influence of elastic 150 i ' : 
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Fig. 4.10 Influence of the 160 
elastic modulus of the (280 MPa, 147 kN) 
confining medium on the 
peak load of rock bolts 1204 J 
z 
p 
& 80+ J 
x 
s 
Q 
a 
v 
= 
[=] 
40} | 
0 , , , , 
0 200 400 600 800 1000 


Young's modulus of the confining medium (MPa) 


4.5.2 Shear Strength of the Bolt/Grout Interface 


To study the influence of the shear strength of the bolt/grout interface on the load 
transfer performance of rock bolts, three analytical pull-out tests were conducted. In 
this case, the shear strength of the bolt/grout interface was varied from 2 to 6 MPa. 
As for the other parameters, they were kept constant. 

The mechanical properties for the rock bolt and the confining medium are tabu- 
lated in Table 4.7. As for the input parameters for the bond-slip model, they are 
tabulated in Table 4.8. 

A plot of the input bond-slip model is shown in Fig. 4.11. Based on the input 
bond-slip model, the calculated load transfer performance of rock bolts is shown in 
Fig. 4.12. 

The results show that the shear strength of the bolt/grout interface played a signif- 
icant role in determining the load transfer performance of rock bolts. Specifically, 
when the shear strength of the bolt/grout interface was 2 MPa, the rock bolt had a 
maximum pull-out load of 146 kN when the pull-out displacement reached 3.1 mm. 


Table 4.7 Mechanical properties of the bolt and the confining medium 
rp (mm) Ep (GPa) L (m) Em (GPa) Am (m°) 
10 196 1.5 1 1 


Table 4.8 Input parameters 


MP. ô MP ô 
for the bond-slip model ga) (Sp (m) Te MPa e D 
2 1.5 0.5 3.5 
4 
6 
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Compared with that, when the shear strength of the bolt/grout interface was 6 MPa, 
the rock bolt had a maximum pull-out load of 282 kN, increasing by 93.2%. And 
the corresponding pull-out displacement reached 4.5 mm. Apparently, the larger 
the shear strength of the bolt/grout interface, the higher the maximum load transfer 
capacity of rock bolts. 

Furthermore, the results show that the shear strength of the bolt/grout interface 
also had an effect in the initial stiffness of the rock bolting system. With the shear 
strength of the bolt/grout interface increasing, the initial stiffness of the rock bolting 
system also increased, as shown in Fig. 4.12. 
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4.5.3 Residual Shear Strength of the Bolt/Grout Interface 


To study the influence of the residual shear strength of the bolt/grout interface on 
the load transfer performance of rock bolts, another rock bolt pull-out scenario was 
created. In this case, the residual shear strength of the bolt/grout interface was varied 
from 0.5 to 1.5 MPa. As for the other parameters, they were kept constant. 

The mechanical properties of the rock bolt and the confining medium are tabulated 
in Table 4.9. As for the input parameters for the bond-slip model, they are tabulated 
in Table 4.10. 

A plot of the input bond-slip model is shown in Fig. 4.13. As for the load transfer 
performance of rock bolts, they are shown in Fig. 4.14. 

The results show that the residual shear strength of the bolt/grout interface had 
no effect on the initial stiffness of the rock bolting system. For example, in this case, 
when the residual shear strength of the bolt/grout interface increased from 0.5 to 
1.5 MPa, the initial stiffness of the rock bolting system kept constant. 

Nevertheless, the residual shear strength of the bolt/grout interface had an effect on 
the maximum load transfer capacity of rock bolts. For example, in this scenario, when 
the residual shear strength of the bolt/grout interface was 0.5 MPa, the maximum pull- 
out load of rock bolts was 226 kN and the corresponding pull-out displacement was 
3.6 mm. Then, when the residual shear strength of the bolt/grout interface increased 
to 1.5 MPa, the maximum pull-out load of rock bolts increased to 256 KN, increasing 
by 13.3%. And the corresponding pull-out displacement was 4.9 mm. Apparently, 
increasing the residual shear strength of the bolt/grout interface was beneficial for 
improving the maximum load transfer capacity of rock bolts. 

Additionally, compared with the shear strength of the bolt/grout interface, 
increasing the residual shear strength of the bolt/grout interface had less apparent 
effect in improving the maximum load transfer capacity of rock bolts. For example, 
in the above analysis, when the residual shear strength of the bolt/grout interface 
increased by 200%, the maximum pull-out load of rock bolts only increased by 
13.3%. However, when the shear strength of the bolt/grout interface increased by 


Table 4.9 Mechanical properties of the bolt and the confining medium 
rp (mm) Ep (GPa) L (m) Em (GPa) Am (m°) 


Table 4.10 Input parameters Tp (MPa) 


ô MP ô 
for the bond-slip model p mm) TEA) : (mm) 


4 1.5 0.5 3:5 
1 
1.5 
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Fig. 4.13 A plot of the input i i ' ENPA 
bond-slip model when the 4} a L 
residual shear strength of the 7=1 MPa 
bolt/grout interface was t =1.5 MPa 
different = NN 
3 b 
Š 
z \ 
a b 
2 
É A 
S 2l S 
z N 
v b T 
= 
N 
1 L 
0 , , , , , 
0 1 2 3 4 5 
Shear slippage (mm) 


Fig. 4.14 Load transfer ' ' ' ' ' 
performance of rock bolts 250} J 
when the residual shear 
strength of the bolt/grout 
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200%, the maximum pull-out load of rock bolts largely increased by 93.2%. There- 
fore, the shear strength of the bolt/grout interface had more significant effect in 
deciding the maximum load transfer capacity of rock bolts. 

The contribution of this paper is that it concentrated on fully grouted rock bolts 
with long embedment length grouted in the surrounding confining medium. Specifi- 
cally, based on the previous work [38, 47], this paper acquired the specific analytical 
solution for this scenario. Furthermore, based on the acquired specific analytical solu- 
tion, it further studied the influence of the elastic modulus of the confining medium 
on the load transfer performance of rock bolts. It revealed that although the peak 
load of rock bolts increased with the elastic modulus of the confining medium, there 
was a critical limit. Once the elastic modulus of the confining medium was beyond 
that critical limit, it had marginal effect on the peak load of rock bolts. 


References 63 


4.6 Conclusions 


A tri-linear model was adopted to depict the bond-slip behaviour of the bolt/grout 
interface. It was assumed that the shear stress of the bolt/grout interface increased 
linearly to the shear strength of the bolt/grout interface with the shear slippage. Then, 
the shear stress decreased linearly to the residual shear strength of the bolt/grout 
interface with the shear slippage. After that, the shear stress was kept constant. 

The pull-out stages of rock bolts were analysed. During the pull-out process, 
the bolt/grout interface underwent elastic, elastic-softening, elastic-softening- 
debonding, softening-debonding and debonding stages. 

Experimental pull-out tests were used to validate the credibility of this analytical 
model. The results show that there was a good correlation between experimental and 
analytical results. 

A parametric study was conducted to evaluate the elastic modulus of the confining 
medium, the shear strength of the bolt/grout interface and the residual shear strength 
of the bolt/grout interface on the load transfer performance of rock bolts. The results 
show that increasing the elastic modulus of the confining medium was beneficial for 
improving the load transfer performance of rock bolts. However, once the elastic 
modulus of the confining medium was beyond a critical limit, it had a marginal 
effect on the peak load of rock bolts. Additionally, increasing the shear strength of 
the bolt/grout interface and the residual shear strength of the bolt/grout interface led 
to rising of the peak load of rock bolts. However, compared with the residual shear 
strength of the bolt/grout interface, increasing the shear strength of the bolt/grout 
interface had more apparent effect in improving the load transfer capacity of rock 
bolts. 
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Chapter 5 A) 
Analytical Studying the Confining get 
Medium Diameter Impact 

on Load-Carrying Capacity of Rock Bolts 


5.1 Introduction 


Rock bolts are tendons or cables that are commonly used in mining engineering 
and civil engineering to reinforce excavated rock masses or soils [1, 2]. Numerous 
in situ tests proved that rock bolts are efficient in improving the internal strength 
of the surrounding confining medium, namely rock masses and soils [3]. Therefore, 
rock bolts are becoming more significant in maintaining the stability of underground 
roadways, chambers, tunnels, surface slopes [4] and rooms in the room and pillar 
mining [5]. In the fully grouted rock bolting system, the bolt can deform consistently 
with the confining medium [6, 7]. Then, the bolt tendon can interact with the confining 
medium to transfer the force [8, 9]. 

To investigate the load-carrying capacity of rock bolts, laboratory pull-out tests 
were the most widely preferred method [10]. Specifically, a natural or artificial rock 
block can be prepared as the confining medium. Then, the bolt was installed in the 
drilled hole in the rock block. Cement-based or resin grout was poured into the 
drilled hole to bond the bolt and rock block. After the grout in the bolting hole 
fully cured, the bolt was pulled out from the rock block. During the experiment, 
load cell and displacement transducers were adopted to monitor the pull-out force 
and displacement. With the recorded data, the force—displacement relationship of 
rock bolts can be plotted, representing the performance of rock bolts. For example, 
Skrzypkowski et al. [11] considered that the load-carrying capacity of segmentally 
installed rock bolts can be influenced by the bolting hole diameter. To investigate this 
impact, they conducted pull-out tests on rock bolts with a diameter of 24 mm. Four 
different bolting hole diameters were used, ranging from 28 to 37 mm. The results 
showed that at the certain installing length of 100 mm and 200 mm, the maximum 
load-carrying capacity of rock bolts varied with the bolting hole diameter. Moreover, 
it was concluded that the optimum load-carrying capacity appeared when the bolting 
hole diameter was 32 mm. 
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The previous research work showed that in the laboratory tests, the geometry of 
the rock block was different. Usually, researchers would like to prefer cylindrical 
or cubic rock blocks. However, it was more common to use cylindrical rock blocks 
[12, 13]. It is interesting to see that although researchers used this method to study 
the performance of rock bolts, the diameter of rock blocks was largely different. For 
laboratory tests, the diameter of rock blocks usually ranged from 100 to 500 mm 
[14]. The previous researchers conducted a number of tests to reveal the loading 
performance of rock bolts when the condition of the bolt, grout and rock masses 
varied [15-17]. However, less attention was paid to evaluate the confining medium 
diameter effect on load-carrying capacity of rock bolts. 

Rajaie [18] was the pioneer to study the confining medium diameter effect. More- 
over, he used concrete materials to cast cylindrical samples. Five different sample 
diameters were used: 100 mm, 150 mm, 200 mm, 250 mm and 300 mm. The results 
showed that when the confining medium diameter was small, the peak force of rock 
bolts ascended rapidly with the confining medium diameter ascending. However, it 
almost stopped ascending once the diameter was beyond 200 mm. This indicated 
that there was a critical influence diameter. Once the rock block diameter was larger 
than this critical influence diameter, further ascending the rock block diameter had 
marginal effect on the performance of rock bolts. Chen et al. [19] continued this 
experimental work, finding that the boundary condition had an effect in determining 
the critical influence diameter. 

Although numerous analytical models were proposed to study the performance of 
rock bolts [20-23], little work has been focused on the confining medium diameter 
effect. Therefore, this paper conducted an analytical study to investigate the influence 
of the confining medium diameter on the load-carrying capacity of rock bolts. First, 
the analytical model concept was illustrated. Then, the confining medium diameter 
effect was studied with this analytical model, and the results were compared with 
experimental results. Last, the influence of the confining medium modulus on the 
critical influence diameter was analysed. 


5.2 Analytical Modelling Approach 


The authors of this paper previously deduced an analytical model to study the loading 
performance of rock bolts [24]. Specifically, after the rock bolt was loaded, shearing 
stress occurred at the bonding face between the tendon and grout. Moreover, shearing 
stress occurred at the bonding face between the grout and confining medium, as shown 
in Fig. 5.1. 

Then, a tri-linear formula was used to illustrate the bonding and debonding 
behaviour of the bonding face, as shown in Eqs. (5.1a-5.1c) [25]: 


t = 23(0 <6 < ô) (5.1a) 
dp 
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Fig. 5.1 State of the bolting z ares eee ee h. ro 
system after the bolt was Pie > 5 confining medium‘: t-t 
loaded A oe ae : 2o gr Bua. 
rock bolt 
~ -= 4 
Tpôr — T,Ô Typ — T, 
i per r“p P TS (Sp <8 <4) (5.1b) 
ô= p 5 — dp 
T= 1,(6 > ôr) (5.1c) 


Then, through analysing the equilibrium relationship between the tensile stress 
in the tendon and shearing stress at the bonding face, the governing equation for the 
bolting system was obtained [24]: 


d2 
sw -xr =0 (5.2) 
4/1 Ap 
v= 
Dy (= t =) G al 


The whole pull-out process of rock bolts was divided into five different phases, 
including the elastic phase (I), the elastic-softening phase (II), the elastic-softening- 
debonding phase (III), the softening-debonding phase (IV) and the debonding phase 
(V) [26]. Through substituting the tri-linear formula and boundary condition in each 
phase into Eq. (5.2), the relationship between the pull-out force and displacement 
can be obtained. 

Specifically, in Phase I, the pull-out force and displacement can be computed as 
follows: 


Dyt, tanh(A, L 
F= Pi eae) (5.4) 
A1Sp 
T, 
i = P52 (5.5) 
Sp 


In Phase II, the pull-out force and displacement can be computed as follows: 


Às us Tp Sp — Trs 
up = = tanh(Ay(L — as)) sin(Azas) — SP cos(azas) + P+? ć (5.6) 
ha A> Tp — T, 


p r 
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1 1 
F= zDot( = tanh(A,(L — a,)) cos(àz2as) + ra sin oa.) ) (5.7) 
1 2 


Tp — Tr 


i 5.8 
3 (5.8) 


x= 


In Phase III, the pull-out force and displacement can be computed as follows: 


ait, le Tak ; 
Up = Sr + + | — tanh (à; (L — ag — a,)) cos(A2a,) + sin(A2as) )da 
2 


2 Ay 
(5.9) 


F= nDs( = tanh(A1(L — aa — ag) cos(àzas) + 2 sin(àza;) + roa) (5.10) 
1 2. 


In Phase IV, the pull-out force and displacement can be computed as follows: 


tan(Ao(L — 
TESE tna S g a 2) (5.11) 
2 iS 
tan(Ao(L — 
F= 1 pyr ( eE Ẹ a) (5.12) 
2, 


Last, in Phase V, the pull-out force and displacement can be computed as follows: 


XGL? 
F=xDyt,( L +8 + h (5.13) 


Equations from (5.4) to (5.13) contributed to the force—displacement relationship 
of rock bolts. Moreover, a good agreement was obtained between the model and 
physical test data [24]. 


5.3 Modelling Process and Results 


To investigate the effect of the confining medium diameter on the load-carrying 
capacity of rock bolts, analytical tests were performed. Specifically, it was assumed 
that pull-out tests were conducted on deformed rock bolts with the diameter of 25 mm 
and elastic modulus of 200 GPa. It was admitted that the deformed rock bolts had 
ribs on the tendon surface. Moreover, the mechanical parameters of ribs, namely rib 
height and rib spacing, may be dependent on the bolt type [27, 28]. However, for 
this paper, at the current stage, there was no specific evaluation on the rib height and 
rib spacing. It was assumed that after the tendon was loaded, the deformed rock bolt 
induced shearing stress at the bonding face between the tendon and grout. As for the 
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Table 5.1 Bonding-slipping MP. MP. 
data for the bonding face one sm Ta MPa) s3 mm) 
between the tendon and grout 2 1 1 3 


influence of the rib height and rib spacing on the shearing stress at the bonding face 
between the tendon and grout, further work will be continued. 

The bolt was installed in the bolting hole with a diameter of 40 mm and grouted 
length of 2 m. For the grout, it was assumed that it had an unconfined compressive 
strength of 50 MPa. For the confining medium, it was assumed that it had an elastic 
modulus of 20 GPa and unconfined compressive strength of 30 MPa. For the bonding 
face between the tendon and grout, the bonding-slipping data shown in Table 5.1 was 
used. 

Table 5.1 summarised the mechanical properties of the bonding face between 
the tendon and grout: shearing strength, slipping where shearing strength reached, 
residual shearing strength and slipping where residual shearing strength reached. 
They can be computed with the following equation [29]: 


P 
T= — (5.14) 
T D,L 

Specifically, the previous research revealed that if the grouted length of the rock 
bolt was short enough, the shearing stress at the bonding face can be treated uniform 
[14]. Therefore, after conducting pull-out tests on rock bolts with short grouted 
length, the maximum pull-out force and residual pull-out force can be substituted 
into Eq. (5.14) to compute the shearing strength and residual shearing strength. Then, 
the slipping where shearing strength reached equalled the pull-out displacement 
where the maximum pull-out force occurred. The slipping where residual shearing 
strength reached equalled the pull-out displacement where the residual pull-out force 
occurred. 

In this paper, for the diameter of the confining medium, four different values were 
used: 100 mm, 150 mm, 200 mm and 250 mm. The configuration of this analytical 
pull-out test is shown in Fig. 5.2. This analytical pull-out process was generally 
consistent with the traditional pull-out experiment of rock bolts [30], where the 
tendon was pulled with a loading rate of 1 mm/min and the load cell together with 
the linear variable differential transformer were used to measure the pull-out force 
and displacement, respectively. 

Then, the analytical pull-out tests were performed, and the results are shown in 
Fig. 5.3. Apparently, the confining medium diameter had an effect in determining the 
performance of rock bolts. Specifically, there was a positive relationship between the 
load-carrying capacity of rock bolts and the confining medium diameter. Moreover, 
the ascending rate of the peak force had a tendency to decline. For example, when 
the confining medium diameter ranged from 200 to 250 mm, the peak force of rock 
bolts was quite close. 
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Fig. 5.2 Configuration of 
the rock bolt loading system 
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Fig. 5.3 Rock bolt performance with different confining medium diameters: a bolting system 
geometry; b force—displacement relationship of rock bolts 


To further investigate the influence of the confining medium diameter on the 
performance of rock bolts, a series of pull-out tests was conducted. In this computa- 
tion, the confining medium diameter varied in the domain [50 mm, 1000 mm]. Then, 
the variation trend of the bolt peak force, when the confining medium diameter was 
varying, was computed, as shown in Fig. 5.4. The peak force ascended significantly 
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when the confining medium diameter varied from 50 to 300 mm. However, when the 
confining medium diameter was larger than 300 mm, the ascending rate of the peak 
force became gentle. Moreover, when the diameter was beyond 400 mm, the peak 
force became almost constant. 

To confirm the accuracy of the analytical modelling results, the experimental tests 
performed by Rajaie [18] were used as a comparison, as shown in Fig. 5.5. It can be 
seen that the peak force variation trend in analytical modelling generally agreed well 
with the experimental test trend. In experimental tests, the peak force rose rapidly 
with the confining medium diameter when it was smaller than 200 mm. After it was 
beyond 250 mm, the peak force became almost unchangeable. This variation trend 
was successfully reflected in analytical modelling. 
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Fig. 5.6 Peak force 
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In the above analysis, the elastic modulus of the confining medium was set as 
20 GPa. However, in experimental tests, the material with different modulus can be 
used to cast the confining medium. Therefore, further pull-out tests were performed 
with different confining medium moduli. Specifically, three different moduli were 
used, 1 GPa, 5 GPa and 20 GPa. The peak force variation trend results are shown in 
Fig. 5.6. Apparently, the peak force variation trend was consistent. This indicated that 
although the confining medium modulus was different, there was a critical influence 
diameter. The confining medium had little effect on the peak force once its diameter 
was over than the critical influence diameter. It was also found that when the confining 
medium modulus was different, the critical influence diameter was also different. For 
example, when the confining medium modulus was 20 GPa, the peak force became 
almost constant when the confining medium diameter was larger than 400 mm. 
Nevertheless, this value was changed when the confining medium with a modulus of 
1 GPa was used. For it, at a diameter of 400 mm, the peak force was still ascending 
rapidly. 

To study the influence of the confining medium modulus on the critical influence 
diameter, further pull-out tests were performed. Specifically, the confining medium 
modulus was ranged from | to 50 GPa. With a certain confining medium modulus, the 
peak force variation trend was first computed. Then, the critical influence diameter 
for that confining medium modulus was recorded. The critical influence diameter 
was determined with the following method. 

First, a relative difference percentage was defined, as shown in Eq. (5.15). 


A= Fmax2 — Fmaxı (5.15) 
Fmaxı 


Where A: the relative difference percentage. F maxı: peak force when the confining 
medium diameter was a. F max2: peak force when the confining medium diameter was 
(a+ Aa). 
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In this paper, the confining medium diameter ascended with a small interval of 
10 mm. Then, when the relative difference percentage was smaller than 0.01%, it 
was assumed that there was no apparent ascending of the peak force. Therefore, 
the corresponding confining medium diameter was regarded as the critical influence 
diameter. 

Following this method, the critical influence diameter for different confining 
medium moduli was acquired, as shown in Fig. 5.7. It can be seen that the confining 
medium modulus had a significant effect in determining the critical influence diam- 
eter. With the confining medium modulus ascending, the critical influence diameter 
declined nonlinearly. This indicated that for the same rock bolt, if the confining 
medium properties were different, the corresponding critical influence diameter may 
be different. Furthermore, the larger the confining medium modulus, the smaller the 
critical influence diameter. 


5.4 Conclusions 


Laboratory pull-out tests were commonly used to study the performance of rock 
bolts. In pull-out tests, natural or artificial cylindrical confining medium was usually 
used to confine rock bolts. This paper conducted an analytical investigation to study 
the influence of the confining medium diameter on the performance of rock bolts. 
In the analytical modelling, the load-carrying capacity of rock bolts with different 
confining medium diameters was compared. The results showed that the confining 
medium diameter had an apparent effect in determining the peak force of rock bolts. 
With the confining medium diameter ascending, the peak force ascended. However, 
the ascending rate of the peak force gradually declined. The peak force variation 
trend was compared with experimental results. And there was a good agreement 
between modelling trend and experimental trend. 
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The rock bolt peak force with different confining medium moduli was studied. 


It was found that although the confining medium modulus was different, the peak 
force variation trend was consistent. For each medium modulus, there was a certain 
critical influence diameter. Beyond that critical influence diameter, further ascending 
the confining medium diameter had marginal effect in improving the peak force of 
rock bolts. 


Moreover, the influence of the confining medium modulus on the critical influence 


diameter was analysed. The results showed that the higher the confining medium 
modulus, the smaller the critical influence diameter. 
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Chapter 6 A) 
Studying the Bond Performance of Fully crest 
Grouted Rock Bolts Based 

on the Variable Controlling 


6.1 Introduction 


Fully grouted rock bolting is commonly used in roadway reinforcement in mining 
and tunnelling reinforcement [1]. Due to the advantage that the reinforcing resistance 
is significant and it can effectively improve the self-bearing capacity of rock masses, 
rock bolting has been widely applied in rock and soil reinforcement [2]. 

However, in situ practice results showed that failure of the reinforcement system 
still occurred [3]. The previous research summarised that failing of the reinforcement 
system could be bolt rupture or bond failure [4]. It is more common to encounter 
bond failure at the boundary surface between bolts and grout [5]. The reason is that 
the bond stress at the boundary surface between bolts and grout is beyond the bond 
strength [6]. More importantly, it is highly related to the bond capacity of bolts [7]. 

To study the bond mechanism of rock bolts, laboratory tests, analytical modelling, 
numerical simulation and engineering practices were used [8]. For example, Zhang 
et al. [9] conducted in situ tests on gob-side roadways to monitor the loading state of 
bolts under the thick coal seam mining condition. Skrzypkowski et al. [10] performed 
laboratory pull-out tests on bolts with different diameters to compare the performance 
of bolts. Chen et al. [11] used the experimental approach to study the anchorage 
performance of bolts installed in a bedded rock mass. Chen and Li [12] adopted the 
modified elements to study the interaction between bolts and in situ rock masses. 

Amongst those methods, analytical modelling has the advantage that it can effec- 
tively reveal the bond course between bolts and rock masses subjected to loading 
[13]. Therefore, this paper studied the bond capacity of bolts based on analytical 
modelling. 

In analytical modelling, former investigators performed numerous investigations 
regarding the bond capacity of bolts [14]. Wang et al. [15] analysed the coupling effect 
of bolts installed in rock masses. Zou and Zhang [16] used a nonlinear equation to 
illustrate the debonding behaviour of the boundary surface between bolts and grout. 
Ren et al. [17] studied the shear failure course of the boundary surface between bolts 
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and grout. Zheng et al. [18] studied the reinforcing system composed of bolts and 
rock masses. Ma et al. [19] established the distributing link of the bond stress at the 
boundary surface between bolts and grout. He et al. [20] deduced the interaction link 
between the tensile stress in bolts and the bond stress after bolts which were loaded. 
Saleem [21] developed an analytical model to study the bond performance of bolts 
under the impact loading condition. The superiority of their model was that it could 
be used to predict the bond performance of different bolt types. Lv et al. [22] deduced 
the load-carrying force of cables based on an elastoplastic model. Wang et al. [23] 
developed a three-dimensional model to study the bond performance of fully grouted 
bolts based on a closed tri-linear relationship. Their results showed that there was 
a critical grouting length. Zhang et al. [24] considered the mutual interaction link 
between bolts and grout. Saleem and Hosoda [25] employed the Latin Hypercube 
sampling approach to analyse the test data regarding several fully grouted bolts. Their 
findings showed that grouting length had more effect on the load-carrying force of 
bolts compared with bolt diameter. 

Those mentioned analytical modelling research successively revealed the coupling 
link between bolts and rock masses. However, much less work has been performed 
to quantitatively study the influence of parameters on the bond performance of fully 
grouted bolts. Moreover, little work was conducted to quantitatively examine the 
bond stress distributing at the boundary surface between bolts and grout. 

To overcome those shortcomings, this paper used the new analytical model devel- 
oped by the authors to quantitatively examine the bond performance of fully grouted 
rock bolts [26]. The originality of this paper is that it used the newly developed 
analytical model to quantitatively examine the bond performance of fully grouted 
bolts. Moreover, relying on this analytical model, the bond stress distributing at the 
boundary surface between bolts and grout was quantitatively studied. 

In this paper, first, the calculating principle of this analytical model was illustrated. 
Then, this analytical model was used to calculate and analyse the load-carrying force 
of bolts. The effect of four different parameters on the load-carrying force of bolts and 
the bond stress distributing state was studied. This study is conducive to enriching 
the base of knowledge. 


6.2 Calculating Principal of the Constitutive Equation 


The authors of this paper previously conducted a theoretical analysis on the load- 
carrying force of bolts. During modelling, the shear course of the boundary surface 
was simulated with a tri-linear equation (Fig. 6.1) [26]. 


t=kô +r (6.1) 


Where k: slope of the tri-linear equation. 
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Fig. 6.1 Tri-linear equation 


Bond stress (MPa) 


Bond slipping (mm) 


The tri-linear equation was revised from the elastic perfectly plastic model and the 
bi-linear model. However, the shortcoming of the elastic perfectly plastic model was 
that it neglected the post-failure performance of the boundary surface. As for the bi- 
linear model, its shortcoming was that it did not consider the residual shear strength 
of the boundary surface [27]. Therefore, through combining those two models, the 
tri-linear equation was used. 

The tri-linear equation assumed that the bond stress of the boundary surface grew 
linearly with the bond slipping. When the bond slipping grew to a specific value, the 
bond stress reached the bond strength. After that, with the continuous growth of the 
bond slipping, the bond stress at the boundary surface decreased to the residual shear 
strength [28]. Last, with the further growth of the bond slipping, the bond stress at 
the boundary surface remained constant, equalling the residual shear strength. 

The credibility of the tri-linear equation was confirmed with the previous exper- 
imental research. For example, Benmokrane et al. [29] conducted pull-out tests on 
fully grouted bolts to study the bond behaviour of the boundary surface. Their results 
showed that after the peak, the shear stress at the boundary surface dropped gradu- 
ally. Moreover, with the slipping further growing, the shear stress was likely to be 
constant. Based on this phenomenon, this paper used the tri-linear equation to depict 
the bond behaviour of the boundary surface. 

Then, this tri-linear equation was incorporated into the anchorage body, and a 
governing equation was developed, as calculated with Eq. (6.2). 


d8? (x) 
dx2 


— r(x) =0 (6.2) 


In this equation, the coefficient of à can be expressed with Eq. (6.3): 
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a Cee (6.3) 
~~ \DbEp  4EmAm 


With the loading of bolts, the boundary surface successively experienced 
five stages, including the elastic stage, the elastic-weakening stage, the elastic- 
weakening-friction stage, the weakening-friction stage and the friction stage [30]. 
Then, analysis was conducted on the tensile stress in bolts and the bond stress on 
the boundary surface during those five stages. Therefore, the equations that describe 
the tensile stress distributing of bolts and the bond stress spreading at the boundary 
surface were obtained. The loading terminal of bolts was regarded as the research 
objective. Then, the link between the load-carrying force and loading displacement 
of bolts was obtained. 

The credibility of this model was validated with experimental tests. The test 
conducted by Bai et al. [31] was used as an example. It showed a good match 
between the test and modelling [26]. 


6.3 Parameter Study 


The input parameters of this model included bolt diameter, bolt modulus, grouting 
length, modulus of rock masses, the cross-sectional area of rock masses and properties 
of the boundary surface. Then, this paper studied the effect of relevant parameters 
on the load-carrying force of bolts. 


6.3.1 Bolt Diameter 


In this section, it was assumed that the bolt diameter was 20 mm, and the modulus 
of bolts was 200 GPa. The grouting length was 2 m. As for the rock masses, the 
modulus was 15 GPa, and the cross-sectional area was 0.5 m2. Those parameters 
were shown in Table 6.1. 

More importantly, the properties of the boundary surface determined the bond 
capacity of bolts. Then, they were shown in Table 6.2. 


Table 6.1 Properties of the anchorage body 
D (mm) E; (GPa) L (m) E> (GPa) A (m?) 
20 200 2 15 0.5 


Table 6.2 Properties of the 


tı (MPa) ô1 (mm) t2 (MPa) ô2 (mm) 
boundary surface 


3 2 1.5 4 
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Fig. 6.2 Effect on bolt 
diameter on the 
load-carrying force of bolts 300+ 


200; 


Force (KN) 
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To study the effect of bolt diameter on the load-carrying force, besides the bolt 
diameter used in Table 6.1, another two bolt diameters were used: 15 mm and 25 mm. 
During the pull-out course, the load-carrying force and loading displacement at the 
loading terminal were recorded. Then, the bond capacity under the loading condi- 
tion was obtained (Fig. 6.2). It shows that bolt diameter significantly affected the 
load-carrying force of bolts. When three different bolt diameters were used, the 
maximum force was 194 KN, 269 kN and 349 KN. With the growth of bolt diameter, 
the maximum force significantly grew. This finding was consistent with the experi- 
mental test results obtained by Kilik et al. [32]. Meantime, it shows that the stiffness 
of the anchorage body grew. When the bolt diameter was 15 mm, the initial stiffness 
of the anchorage body was 50 kN/mm. However, when the bolt diameter grew to 
25 mm, the stiffness of the anchorage body grew to 104 kN/mm. 

Also, it shows that bolt diameter had specific effect on the post-failure performance 
of the reinforcement system. For instance, when bolt diameter was small, such as 
15 mm, after bolts reached the maximum force, the loading displacement showed 
an apparent snapback phenomenon. This was because the boundary surface entered 
the weakening-friction stage. Moreover, it was directly resulted by the continuous 
growth of the friction length. This phenomenon also occurred in the previous research 
[33]. By contrast, when bolt diameter was relatively larger, such as 25 mm, after bolts 
reached the maximum force, this snapback phenomenon was not apparent. 

Additionally, the relationship between the maximum force of bolts and A was 
analysed. With bolt diameter growing, the corresponding à dropped from 3.66e—5 
to 2.85e—5 N/m (Fig. 6.3). As for the maximum force of bolts, it grew directly. 
Therefore, a negative relationship between the maximum force of bolts and À 
occurred. 

This paper mainly focused on the maximum force of bolts. Therefore, when those 
bolts with three different diameters reached the maximum force, the bond stress 
distribution of the boundary surface was exported (Fig. 6.4). It shows that when 
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bolts with those three different diameters reached the maximum force, the bond 
stress at the boundary surface was in the elastic-weakening-friction stage. 

Meantime, it shows that for bolts with a smaller diameter, when the maximum 
force was reached, the length of the friction section was larger. For instance, when the 
bolt diameter was 15 mm, the corresponding friction length was 828 mm. However, 
when the bolt diameter was 25 mm, the friction length dropped to 493 mm. This 
was related to the reinforcement stiffness of the anchorage body. When the bolt 
diameter was small, the reinforcement stiffness was also slight. Therefore, to reach 
the maximum force, a relatively large loading displacement was needed. However, 
under the effect of the large loading displacement, the boundary surface experienced 
rather large bond slipping. Then, the friction length grew. 


Fig. 6.4 Effect of bolt 
diameter on the bond stress 
distributing on the boundary 
surface when peak load 
occurred 


Bond stress (MPa) 


Grouting length (m) 
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6.3.2 Elastic Modulus of Bolts 


In the in situ condition, rock bolts with different modulus may be used. For instance, 
in the permanent reinforcement, steel bolts with a modulus of 200 GPa were always 
used. However, in the roadway along the working face side, the fibre-reinforced 
polymer bolts with modulus of 50 GPa may be used. 

Therefore, the modulus of bolts was regarded as the research objective. The 
purpose was to analyse its effect on the load-carrying force of bolts. In this section, 
a calculation was conducted on bolts with a diameter of 20 mm. To study the effect 
of the modulus of bolts, three pull-out analysis was conducted. The corresponding 
modulus of bolts was 50 GPa, 100 GPa and 200 GPa. The other input parameters 
were the same as the values used in Sect. 6.3.1. The comparing results of those three 
pull-out tests are shown in Fig. 6.5. It shows that the elastic modulus of bolts signifi- 
cantly affected the pre-peak and post-peak performance of fully grouted bolts. First, 
with the modulus of bolts growing from 50 to 200 GPa, the maximum force grew 
from 229 to 269 KN, growing by 17%. This indicated that although increasing the 
modulus of bolts was conducive to improving the bond performance of bolts, this 
impact was not significant. Secondly, the modulus of bolts had a significant impact 
on the initial stiffness of the anchorage body. Specifically, when the modulus of bolts 
was 50 GPa, the corresponding initial stiffness of the anchorage body was 38 kN/m. 
When the modulus of bolts grew to 200 GPa, the corresponding initial stiffness of the 
anchorage body was 76 kN/m, growing by 100%. Therefore, the larger the modulus 
of bolts, the higher the stiffness of the anchorage body. 

Thirdly, the modulus of bolts significantly influenced the loading displacement 
where the maximum force of bolts occurred. For example, when bolts with a modulus 
of 50 GPa were tested, the loading displacement at peak force was 18.1 mm. When the 
modulus of bolts grew to 100 GPa, the loading displacement at peak force dropped to 


Fig. 6.5 Effect of the 300 
modulus of bolts on the 
load-carrying force 


200; 
Z 
Š 
g 
2 
© 
ey 
100; 
—— E=50 GPa 
——— E=100 GPa 
—— E=200 GPa 
0 1 L T 
0 5 10 15 20 


Displacement (mm) 


88 6 Studying the Bond Performance of Fully Grouted Rock Bolts Based ... 


i : ‘ 8 280 
Fig. 6.6 Relationship P 
between the maximum force —— Force 
of bolts, à and elastic 7 1268 
modulus of bolts 
Bg 1256 & 
3 4 
A mi 
S 2 
x 5 244 £ 
“= 
4 232 
3 : : 220 
50 100 150 200 


Elastic modulus of bars (GPa) 


10.6 mm. When bolts with a modulus of 200 GPa were tested, the loading displace- 
ment dropped to 6.49 mm. Therefore, the loading displacement at peak force signif- 
icantly dropped with the modulus of bolts. Last, the modulus of bolts influenced the 
post-peak performance of bolts. When the modulus of bolts was 50 GPa, after the 
peak load, there was a gradual dropping of the load-carrying force. By contrast, when 
the modulus of bolts grew to 200 GPa, after the peak force, the load-carrying force 
dropped dramatically. 

Also, the relationship between the maximum force of bolts and à was analysed. 
In this case, à dropped from 6.33e—5 to 3.18e—5 N/m (Fig. 6.6). By contrast, the 
maximum force of bolts grew. Therefore, for this case, the smaller the 4, the larger 
the maximum force of bolts. 

When bolts reached the maximum force, the bond stress distributing at the 
boundary surface is exported (Fig. 6.7). The modulus of bolts had a significant 
effect on the bond stress distributing at the boundary surface. When the modulus of 
bolts was small, the friction length of the boundary surface was relatively large. For 
instance, when the modulus of bolts was 50 GPa, the friction length of the boundary 
surface was 1322 mm. However, when the modulus of bolts was 200 GPa, the friction 
length of the boundary surface was 648 mm. This was because when the modulus 
of bolts was small, the stiffness of the reinforcement system was relatively small. 
Therefore, after a rather large displacement, bolts finally reached the maximum force. 
Moreover, when a large loading displacement generated, the boundary surface expe- 
rienced relatively large bond slipping. Therefore, the friction length was relatively 
large. 


6.3.3 Grouting Length of Bolts 


In engineering practices, the grouting length can be significantly different [34]. There- 
fore, in this section, the effect of grouting length on the load-carrying force of bolts 
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Fig. 6.7 Effect of the 3 
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surface when peak force 
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was studied. The calculation was conducted on bolts with a diameter of 20 mm 
and modulus of 200 GPa. During the calculating course, the grouting length ranged 
between 2 and 3 m. The other input parameters were the same as the values used in 
Sect. 6.3.1. Then, under those three different grouting length conditions, the bond 
capacity of bolts was acquired (Fig. 6.8). It shows that growth of the grouting length 
did not affect the initial reinforcement stiffness of bolts. Under those three different 
grouting length conditions, in the elastic stage, the initial stiffness was almost same. 
However, the grouting length had a significant effect on the maximum force of bolts. 
With the growth of grouting length, the maximum force of bolts significantly grew. 
When the grouting length grew from 2 to 3 m, the maximum force of bolts grew 
from 269 to 364 kN. 


Fig. 6.8 Effect of the 400 
grouting length on the 
load-carrying force of bolts 
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Attention is paid that this analytical model can be used to calculate the ideal 
grouting length for bolts. For example, in this case, it is assumed that the tensile 
strength of bolts was 1200 MPa. Then, this tensile strength was substituted into 
Eq. (6.4) to calculate the rupture force of bolts. 


or D? 
F, = —— (6.4) 
4 
Where F,: rupture force of bolts. 

The calculating result showed that the rupture force of bolts was 377 KN. Then, this 
rupture force can be substituted into the present analytical model to back-calculate 
the ideal grouting of bolts. For this case, results showed that the ideal grouting length 
of bolts was around 3.5 m, as shown in Fig. 6.9. 

When bolts reached the maximum force, the bond stress at the boundary surface 
was exported (Figs. 6.10, 6.11 and 6.12). Being different from other parameters, 
the grouting length had a marginal effect on the bond stress distributing state at the 
boundary surface. Under three different grouting length conditions, the elastic length 
of the boundary surface was around 552 mm. All the weakening length was around 
800 mm. The only difference was the friction length. Under those three grouting 
length conditions, the friction lengths were 648 mm, 1149 mm and 1648 mm. 

Figure 6.11 showed the bond stress distributing at the boundary surface at the 
peak force when the grouting length was 2.5 m. It clearly indicated that when the 
peak force was reached, the bond stress at the boundary surface was in the elastic- 
weakening-friction state. At the free end of bolts, the bond stress was 2.42 MPa. 
With the distance from the free end growing, the bond stress grew nonlinearly. When 
the distance from the free end grew to 0.55 m, the bond stress reached a peak of 
3 MPa. With the distance from the free end further growing, the bond stress dropped 
nonlinearly. When the distance from the free end grew to 1.35 m, the bond stress 
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Fig. 6.10 Bond stress l i 
distributing at the boundary 
surface at the peak force 
when the grouting length 
was 2 m 
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Fig. 6.11 Bond stress 
distributing at the boundary 
surface at the peak force 
when the grouting length 
was 2.5 m 
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dropped to the residual bond strength of 1.5 MPa. After that, the bond stress remained 
constant, levelling at 1.5 MPa. 

Figure 6.12 showed the bond stress distributing at the boundary surface at the 
peak force when the grouting length was 3 m. It showed that although the grouting 
length grew, the bond stress distributing was still in the elastic-weakening-friction 
state when the maximum force was reached. More interestingly, the positions where 
the bond strength and the residual bond strength were reached were identical with 
Fig. 6.11. However, the only difference was that the friction length grew. In this case, 
the friction length grew to 1648 mm, 499 mm longer than the friction length shown 
in Fig. 6.11. 
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Fig. 6.12 Bond stress 3L [— 1 m | 
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6.3.4 Bond Slipping When the Bond Strength Was Reached 


Considering that the properties of the boundary surface affected the bond capacity of 
bolts, the bond slipping at the bond strength was regarded as a factor to conduct the 
analysis. During the studying, a calculation was conducted on a bolt with a diameter 
of 20 mm, modulus of 200 GPa and grouting length of 2 m. To study the effect of the 
bond slipping, it ranged between | and 3 mm. The corresponding tri-linear equation 
is shown in Fig. 6.13. 

During calculating, the other parameters were the same as the values used in 
Sect. 6.3.1. After that, the analytical pull-out test was conducted on bolts, and the 
corresponding load-carrying force was obtained (Fig. 6.14). It shows that with the 
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growth of the bond slipping at the bond strength, the maximum force of bolts 
decreased gradually. Meantime, the stiffness of the reinforcing system gradually 
decreased. However, this bond slipping did not affect the residual load-carrying force 
of bolts. Therefore, under those three different conditions, the residual load-carrying 
force of bolts remained around 188 KN. 

When the maximum force of bolts was reached, the bond stress at the boundary 
surface was exported (Fig. 6.15). With the growth of the bond slipping at the bond 
strength, the elastic length and the friction length of the boundary surface gradually 
grew. For instance, as the bond slipping at the bond strength grew from 1 to 3 mm, 
the elastic length of the boundary surface grew from 300 to 785 mm. Moreover, the 
corresponding friction length of the boundary surface grew from 483 to 814 mm. 
The weakening length dramatically declined from 1218 to 400 mm. 
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6.4 Limitation of the Current Study 


The limitation of the current paper was that it neglected the influence of the filling 
grout in the borehole. It only considered the interaction between fully grouted bolts 
and surrounding rock mass. 


6.5 Recommendation of Further Work 


After the current study, the authors plan to incorporate this analytical model into the 
underground reinforcement scenarios. Then, the bond performance of fully grouted 
bolts reinforced in underground roadways and tunnels can be studied. In this way, the 
reinforcement parameters of fully grouted bolts used in the underground roadways 
and tunnels can be optimised. 


6.6 Conclusions 


This paper conducted an analytical study on the bond performance of fully grouted 
rock bolts based on the variable controlling method. Relied on this approach, the 
effect of bolt diameter, modulus of bolts, grouting length and bond slipping at the 
bond strength on the load-carrying force of bolts were studied. Moreover, the bond 
stress distributing at the boundary surface was studied. This paper is conducive to 
enriching the basement of knowledge. The main conclusions are listed below: 


(1) With bolt diameter growing, the maximum force of bolts significantly grew. In 
this case, when bolt diameter grew from 15 to 25 mm, the maximum force of 
bolts grew from 194 to 349 KN, growing by 80%. Meantime, the stiffness of the 
reinforcement system grew by 108%. When bolts reached the maximum force, 
the bond stress at the boundary surface was in the elastic-weakening-friction 
state. Moreover, the smaller the bolt diameter, the longer the friction length. 

(2) The elastic modulus of bolts significantly affected the load-carrying force of 
bolts. With the growth of bolts’ modulus, the maximum force of bolts gradually 
grew. When the elastic modulus grew from 50 to 200 GPa, the maximum force 
grew from 229 to 269 kN, only growing by 17%. However, the stiffness of 
the reinforcement system grew by 100%. Moreover, the friction length at the 
boundary surface declined. 

(3) The grouting length had a significant effect on the maximum force of bolts. 
With the grouting length growing, the maximum force of bolts apparently grew. 
With the grouting length growing from 2 to 3 m, the maximum force rose from 
269 to 364 kN, growing by 35%. However, the grouting length had almost 
no effect on the stiffness. When bolts reached the maximum force, the bond 
stress distributing of the boundary surface was consistent. With the grouting 
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length growing, the length of the elastic part and the weakening part was equal. 
By contrast, only the friction length grew. With the tensile strength of bolts 
provided, the current analytical model can be used to back-calculate the ideal 
grouting length of bolts. 


(4) The bond slipping at the bond strength affected the load-carrying force of bolts. 
When it grew from 1 to 3 mm, the maximum force dropped from 281 to 258 kN, 
dropping by 8%. As for the stiffness, it dropped by 44%. Moreover, the elastic 
length and the friction length grew. By contrast, the weakening length declined. 
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Chapter 7 A) 
Analytical Modelling Rock Bolts rie 
with a Closed Nonlinear Model 


7.1 Introduction 


Rock bolts have been used in the mining industry for a long time. They are used as 
the reinforcement tendons to keep the stability of openings in underground mining 
and the slope in open-pit mining. Furthermore, currently, rock bolts are commonly 
used not only in the mining engineering but also in civil engineering [1]. 

The rock bolt is a tendon which can be a rod or a cable inserted into a bore- 
hole drilled in rock masses or soils. Traditionally, the material of rock bolts was 
steel. However, engineering practices showed that the steel rock bolts can be easily 
corroded by the moist environment in underground mines. Then, the fibre-reinforced 
polymer (FRP) was used to manufacture rock bolts, which are termed as FRP rock 
bolts. Experimental results proved that the FRP rock bolts can effectively reinforce 
the surrounding rock masses without being corroded [2]. 

When rock bolts are installed, they can be point anchored or fully grouted. And 
the bonding agent can be polyester resin or cementitious grout. This chapter deals 
with fully grouted rock bolts. 

After fully grouted rock bolts are installed in rock masses, the movement of 
surrounding rock mass induces shear deformation of the grout column. Consequently, 
shear stress occurs along the bolt/grout and grout/rock interfaces. Then, the load 
can be transferred between rock bolts and the surrounding rock masses. To study 
the load transfer process of rock bolts, numerous laboratory and in situ tests have 
been conducted. Compared with that, relatively less work has been conducted with 
analytical modelling to study the load transfer behaviour of rock bolts. 

The pioneering work in this area should be attributed to Hawkes and Evans [3], 
who developed an analytical model to study the shear stress distribution along a 
fully grouted steel bolt. However, their results were only applicable to the case that 
the grouted bolt was coupled with the surrounding confining medium. Farmer [4] 
studied the shear stress distribution of the bolt/grout interface with an analytical 
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model. Experimental work was used to validate this analytical model. However, this 
model was not applicable for the debonding behaviour of the bolt/grout interface. 

Ren et al. [5] studied the load transfer behaviour of a fully grouted rock bolt 
with a tri-linear model. The pull-out stages of a rock bolt were divided into five 
stages, including elastic, elastic-softening, elastic-softening-debonding, softening- 
debonding and debonding. The model was successfully validated with laboratory 
and in situ test results. However, when the tri-linear model was used to depict the 
bond-slip behaviour of the bolt/grout interface, pure softening may also occur along 
the bolt/grout interface, as indicated by Blanco Martin et al. [6]. And this was not 
considered in the analytical model proposed by Ren et al. [5]. 

Recently, Liu et al. [7] proposed an analytical model to study the load transmission 
model for a rock bolt subjected to open and sliding joint displacements. Experimental 
work was used to validate this model. The modelling results revealed that for a definite 
joint displacement and joint location, larger anchored angle between the rock bolt 
and the joint face was beneficial for improving the reinforcement effect. 

Ma et al. [8] proposed an analytical model to study the load transfer mechanism 
of fully grouted rock bolts. In their model, a closed nonlinear model was used to 
depict the bond-slip behaviour of the bolt/grout interface [8—10]. In this chapter, the 
authors adopted their model and conducted a study on the load transfer behaviour of 
rock bolts. First, a general illustration of the analytical model was given. Then, an 
experimental test was used to validate this analytical model. Following it, a parametric 
study was conducted with this analytical model. Last, the load distribution along a 
fully grouted rock bolt was studied based on this analytical model. 


7.2 Illustration of the Analytical Model 


For a fully grouted rock bolting system without faceplate, when the rock bolt is 
loaded, shear stress occurs along the bolt/grout and grout/rock interfaces. When the 
shear stress is higher than the interfacial bond strength, bond failure of the interface 
occurs. Experimental tests proved that the major failure mode is the bond failure of 
the bolt/grout interface [11]. Therefore, in this chapter, failure along the bolt/grout 
interface was studied. 

Assuming that the rock bolt is in linear elastic deformation along the axial direction 
and the shear stress of the bolt/grout interface is equal to the axial displacement of 
the rock bolt at the same position, the shear stress of the bolt/grout interface, the 
pull-out load of the rock bolt and the load distribution along the rock bolt can be 
expressed with Eqs. (7.1), (7.2) and (7.3), respectively [8, 9]. 


1(8) = mee Femi - e7) (7.1) 


Where a and b are coefficients. 
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Where xo is a coefficient and can be calculated with Eq. (7.4). 
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Xo = n{ — — : 
0 b 4h, 


7.3 Validation of the Analytical Model 


To further confirm the credibility of this analytical model, an in situ pull-out test was 
used to validate this analytical model. Bastamia et al. [12] conducted a case study in a 
coal mine. Specifically, rock bolts with a diameter of 22 mm were installed in the rock 
mass [12]. The rock bolt has an elastic modulus of 207 GPa and embedment length 
of 1800 mm. During the test, the pull-out load and displacement were recorded. 

In this chapter, to analytically model the pull-out process of the rock bolt, Eq. (7.2) 
was used to calculate the load versus displacement of the rock bolt. It was found 
that when a = 1.5 mm and b = 620.6 mm, there was a close match between the 
experimental pull-out test result and the analytical result, as shown in Fig. 7.1. 

By substituting the calculated coefficients a and b into Eq. (7.1), the shear stress 
versus shear slippage of the bolt/grout interface can be acquired, as shown in Fig. 7.2. 


Fig. 7.1 Comparison of the 
in situ test result conducted 
by Bastamia et al. [12] and 
the analytical modelling 
result 
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Fig. 7.2 Shear stress versus 
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It can be seen that after the rock bolt was loaded, the shear stress of the bolt/grout 
interface increased nonlinearly. When the shear slippage of the bolt/grout interface 
was 1 mm, the shear stress of the bolt/grout interface reached the bond strength of the 
bolt/grout interface, which was 1.11 MPa. After that, the shear stress of the bolt/grout 
interface decreased exponentially. 


7.4 Parametric Study 


After the analytical model was validated, a parametric study was conducted with this 
analytical model to evaluate the influence of different parameters on the load transfer 
behaviour of rock bolts. 


7.4.1 The First Coefficient 


In this analytical model, there are two important coefficients including a and b. The 
parametric study was first conducted to evaluate the influence of the coefficient a on 
the load transfer performance of rock bolts. Specifically, the load transfer of a rock 
bolt with an elastic modulus of 210 GPa and diameter of 20 mm was simulated. The 
coefficient b was kept constant in this calculation and equal to 800 mm. As for the 
coefficient a, three different values were used separately, namely 2, 3 and 4 mm. The 
pull-out load versus displacement relationship is shown in Fig. 7.3. 

Apparently, the coefficient a had an obvious effect on the load transfer capacity 
of rock bolts. The maximum load transfer capacity of the rock bolt was extracted. 
Then, the relationship between the maximum load transfer capacity of the rock bolt 


7.4 Parametric Study 103 


Fig. 7.3 Pull-out load 
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and the coefficient a can be plotted, as shown in Fig. 7.4. It can be seen that with the 
coefficient a increasing, the maximum load transfer capacity of rock bolts increases 
linearly. 

Meanwhile, the shear stress versus shear slippage of the bolt/grout interface can 
be acquired, as shown in Fig. 7.5. 

The results show that although the coefficient a was different, the trend of the bond- 
slip behaviour of the bolt/grout interface was still consistent. However, the coefficient 
a had an apparent influence in deciding the shear slippage when the bond strength 
of the bolt/grout interface occurred. For example, with the coefficient a increasing 
from 2 to 4 mm, when the bond strength of the bolt/grout interface occurred, the 
shear slippage also increased from 1.4 to 2.8 mm. 
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Fig. 7.5 Shear stress versus 2mm 
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Furthermore, with the coefficient a increasing, the bond strength of the bolt/grout 
interface also increased. To further illustrate this, the bond strength of the bolt/grout 
interface when the coefficient a was different was extracted. Then, the relationship 
between the bond strength of the bolt/grout interface and the coefficient a was plotted, 
as shown in Fig. 7.6. 

The results show that when a = 2 mm, the bolt/grout interface had a bond strength 
of 0.82 MPa. Compared with that, when a = 3 mm, increasing by 50%, the bolt/grout 
interface had a bond strength of 1.23 MPa, increasing by 50%. Last, when a = 
4mm, increasing by 100%, the bolt/grout interface had a bond strength of 1.64 MPa, 
increasing by 100%. Apparently, the larger the coefficient a, the higher the bond 
strength of the bolt/grout interface. Additionally, there was a linear relationship 
between the bond strength of the bolt/grout interface and the coefficient a. 


Fig. 7.6 Effect of coefficient 
a on the bond strength of the 
bolt/grout interface 1.54 J 


Bond strength (MPa) 


0.5F 


3 4 


Value of the coefficient a (mm) 


7.4 Parametric Study 105 


7.4.2 The Second Coefficient 


Similarly, the influence of coefficient b on the load transfer performance of a rock 
bolt was studied. The rock bolt had an elastic modulus of 210 GPa and diameter of 
20 mm. In this case, the coefficient a was kept constant and equal to 2 mm. As for the 
coefficient b, three different values were used, namely 500, 600 and 700 mm. The 
pull-out load versus displacement relationship of the rock bolt is shown in Fig. 7.7. 
Apparently, the coefficient b had a significant effect on the load transfer perfor- 
mance of rock bolts. With the coefficient b increasing, the load transfer capacity 
of rock bolts tended to decrease. To clearly show this, the maximum load transfer 
capacity of rock bolts was extracted and plotted with the coefficient b (Fig. 7.8). 
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Fig. 7.9 Shear stress versus b=500mm 
shear slippage curve of the 2L — b=600mm | 
bolt/grout interface when —— b=700mm 
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It shows that when b = 500 mm, the maximum load transfer capacity of the rock 
bolt was 264 kN. Compared with that, when b = 600 mm, increasing by 20%, the 
maximum load transfer capacity of the rock bolt was 220 kN, decreasing by 16.7%. 
When b = 700 mm, increasing by 40%, the maximum load transfer capacity of the 
rock bolt was 188 kN, decreasing by 28.8%. Therefore, increasing the coefficient b 
leaded to nonlinear decreasing of the maximum load transfer capacity of rock bolts. 
As for the influence of the coefficient b on the bond-slip behaviour of the bolt/grout 
interface, it is shown in Fig. 7.9. 

It shows that with the coefficient b increasing, the trend of the bond-slip of the 
bolt/grout interface kept consistent. Furthermore, the coefficient b had no effect 
in deciding the shear slippage where the bond strength of the bolt/grout interface 
occurred. For example, in this case, although the coefficient b was different, the 
shear stress of the bolt/grout interface reached the bond strength at the same shear 
slippage of 1.4 mm. 

However, the coefficient b had an apparent effect on the bond strength of the 
bolt/grout interface. To illustrate this clearly, the bond strength of the bolt/grout 
interface was extracted and plotted with the coefficient b, as shown in Fig. 7.10. 

The results show that increasing the coefficient b leaded to decreasing the bond 
strength of the bolt/grout interface. Specifically, when b = 500 mm, the bond 
strength of the bolt/grout interface was 2.1 MPa. Compared with that, when b = 
600 mm, increasing by 20%, the bond strength of the bolt/grout interface decreased 
to 1.46 MPa, decreasing by 30.5%. When b = 700 mm, increasing by 40%, the 
bond strength of the bolt/grout interface decreased to 1.07 MPa, decreasing by 49%. 
Therefore, there was a negative nonlinear relationship between the bond strength of 
the bolt/grout interface and the coefficient b. 
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Fig. 7.10 Effect of 
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7.4.3 Elastic Modulus of the Rock Bolt 


The influence of the elastic modulus of the rock bolt was studied on the load transfer 
performance of rock bolts. In this case, the rock bolt had a diameter of 20 mm. As 
for the coefficients a and b, they were 2 and 500 mm. As for the elastic modulus of 
the rock bolt, they were ranged from 70 to 210 GPa with an interval of 70 GPa. The 
load versus displacement relationship of the rock bolt is shown in Fig. 7.11. 

The results show that increasing the elastic modulus of the rock bolt was beneficial 
for improving the maximum pull-out load of the rock bolt. Furthermore, with the 
elastic modulus of the rock bolt increasing, the slope of the load—displacement curve 
also increased. This indicated that increasing the elastic modulus of the rock bolt 
was beneficial for improving the stiffness of the rock bolting system. 


Fig. 7.11 Pull-out load 
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Fig. 7.12 Pull-out load 
versus displacement curve of 400; 
the rock bolt when the rock 
bolt diameter was different 
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7.4.4 Rock Bolt Diameter 


The analytical model was also used to evaluate the influence of the rock bolt diameter 
on the load transfer behaviour of rock bolts. Specifically, the elastic modulus of the 
rock bolt was kept constant and equal to 210 GPa. As for the coefficients a and b, 
they were set as 2 mm and 500 mm. The rock bolt diameter was ranged from 15 to 
25 mm with an interval of 5 mm. The pull-out load versus displacement relationship 
of rock bolts is shown in Fig. 7.12. It shows that increasing the rock bolt diameter 
can effectively improve the load transfer capacity of rock bolts. 


7.5 Load Distribution Along a Fully Grouted Rock Bolt 


The load distribution along a fully grouted rock bolt was studied with this analytical 
model. Specifically, the rock bolt had an elastic modulus of 210 GPa and diameter 
of 25 mm. As for the coefficients a and b, they were set as 0.53 mm and 200 mm. 
To simulate the load variation along the grouted rock bolt, an embedment length of 
2000 mm was used. 

The rock bolt was pulled to a displacement of 3.2 mm and reached a maximum 
pull-out load of 273 kN. During the pull-out process, the axial load in the rock bolt 
was calculated when the pull-out load ranged from 50 to 250 kN with an interval of 
50 KN. Specifically, by substituting the pull-out load to Eq. (7.4), xo can be acquired. 
Then by substituting xo to Eq. (7.3), the axial load along the rock bolt from the free 
end to the loaded end can be acquired, as shown in Fig. 7.13. 

The results show that during the pull-out process, the axial load in the rock 
bolt decayed from the loaded end to the free end independent of the pull-out load. 
However, the decaying trend of the axial load was different. Specifically, when the 
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pull-out load was small, for example, 50 KN, the axial load in the rock bolt decayed 
from the loaded end to the free end with a simple exponential form. Furthermore, the 
bending direction of the load distribution curve was always upward. However, when 
the pull-out load was large, for example, 250 KN, only the axial load in the rock bolt 
close to the free end decayed with an exponential form. There were two bending 
directions of the load distribution curve. For the section close to the free end, the 
bending direction of the load distribution curve was upward whilst for the section 
close to the loaded end, the bending direction of the load distribution curve was 
downward. Nevertheless, further experimental work should be conducted to validate 
this finding. 


7.6 Limitation and Recommendation for Future Work 


A limitation of the current study was that only one experimental pull-out test was 
used to validate the credibility of this analytical model. In fact, the previous research 
has proved that the rock bolt performance was influenced by many parameters such as 
the rock bolt profile, the rock mass property, the embedment length and the borehole 
diameter [13-18]. Therefore, in the future work, more experimental pull-out tests 
with different performance will be adopted to further confirm the accuracy of the 
analytical model. 
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7.7 Conclusions 


An analytical model was used to study the load transfer behaviour of fully grouted 
rock bolts. In this analytical model, a closed nonlinear model was used to depict the 
bond-slip behaviour of bolt/grout interface. 

An in situ pull-out test was used to validate this analytical model, showing 
that there was a close match between the experimental result and analytical result. 
Following it, a parametric study was conducted to evaluate the influence of two coef- 
ficients, elastic modulus of the rock bolt and the diameter of the rock bolt on the 
performance of rock bolts. It was found that increasing the coefficient a leaded to 
linear increasing of the maximum pull-out load of rock bolts and the bond strength 
of the bolt/grout interface. On the other hand, increasing the coefficient b leaded 
to nonlinear decreasing of the maximum pull-out load of rock bolts and the bond 
strength of the bolt/grout interface. Improving the elastic modulus of the rock bolt 
and the diameter of the rock bolt were beneficial for improving the load transfer 
capacity of the rock bolting system. 

Additionally, the load distribution along the fully grouted rock bolt was studied 
with this analytical model. It was found that axial load in the rock bolt decayed 
from the loaded end to the free end independent of the pull-out load. However, the 
trend of the load distribution curve along the rock bolt was influenced by the pull-out 
load. When the pull-out load was relatively small, the bending direction of the load 
distribution curve was upward. However, when the pull-out load was relatively large, 
there were two bending directions of the load distribution curve. For the section close 
to the free end, the bending direction of the load distribution curve was upward whilst 
for the section close to the loaded end, the bending direction was downward. 
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